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Abstract 

This  thesis  investigates  a  modification  to  Differential  Phase  Shift  Keyed  (DPSK) 
modulation  to  create  a  Low  Probability  of  Interception/Exploitation  (LPI/LPE)  com¬ 
munications  signal.  A  pseudorandom  timing  offset  is  applied  to  each  symbol  in  the 
communications  stream  to  intentionally  create  intersymbol  interference  (ISI)  that 
hinders  accurate  symbol  estimation  and  bit  sequence  recovery  by  a  non-cooperative 
receiver.  Two  cooperative  receiver  strategies  are  proposed  to  mitigate  the  ISI  due  to 
symbol  timing  offset:  a  modified  minimum  Mean  Square  Error  (MMSE)  equalization 
algorithm  and  a  multiplexed  bank  of  equalizer  filters  determined  by  an  adaptive  Least 
Mean  Square  (LMS)  algorithm.  Both  cooperative  receivers  require  some  knowledge 
of  the  pseudorandom  symbol  timing  dither  to  successfully  demodulate  the  communi¬ 
cations  waveform. 

Numerical  Mat  lab®  simulation  is  used  to  demonstrate  the  bit  error  rate  per¬ 
formance  of  cooperative  receivers  and  notional  non-cooperative  receivers  for  binary, 
4-ary,  and  8-ary  DPSK  waveforms  transmitted  through  a  line-of- sight,  additive  white 
Gaussian  noise  channel.  Simulation  results  suggest  that  proper  selection  of  pulse 
shape  and  probability  distribution  of  symbol  timing  offsets  produces  a  waveform  that 
is  accurately  demodulated  by  the  proposed  cooperative  receivers  and  significantly 
degrades  non-cooperative  receiver  symbol  estimation  accuracy.  In  typical  simula¬ 
tions,  non-cooperative  receivers  required  2-8  dB  more  signal  power  than  cooperative 
receivers  to  achieve  a  bit  error  rate  of  1.0%.  For  nearly  all  reasonable  parameter 
selections,  non-cooperative  receivers  produced  bit  error  rates  in  excess  of  0.1%,  even 
when  signal  power  is  unconstrained. 
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Signal  Processing  Design  of 


Low  Probability  of  Intercept  Waveforms 
Via  Intersymbol  Dither 

I.  Introduction 

Secure  communications  have  proven  to  be  vital  to  successful  military  operations 
across  all  levels  of  warfare  [9].  With  the  advent  of  wireless  technology,  dis¬ 
tributed  communication  has  become  a  military  staple  [17].  The  advantage  of  wireless 
communications  is  not  without  risks.  In  particular,  the  broadcast  nature  of  wireless 
transmissions  makes  this  form  of  communication  highly  susceptible  to  eavesdropping 
by  third  parties  utilizing  non-cooperative  receiver  devices.  This  vulnerability  has  led 
to  an  ever-evolving  challenge  to  develop  and  field  wireless  communications  systems 
that  hinder  or  preclude  detection,  interception,  and  exploitation  by  non-cooperative 
receivers  [12]. 

A  multitude  of  strategies  have  been  employed  to  prevent  the  exploitation  of 
communication  signals.  Perhaps  the  most  common  approach  is  data  encryption  which 
involves  coding  of  the  information  bits  prior  to  transmission.  Successful  encryption 
prevents  a  third  party  from  extracting  information  from  the  intercepted  and  demod¬ 
ulated  communications  waveform  [21:463]. 

Another  way  to  approach  wireless  communications  security  is  to  prevent  the 
adversary  from  intercepting  the  waveform  in  a  meaningful  way.  Spread  spectrum  and 
frequency  hopping  techniques,  for  example,  enhance  wireless  communications  security 
by  making  the  signal  difficult  to  detect  and  intercept  [21:95].  This  thesis  investigates  a 
strategy  that  manipulates  the  transmitted  waveform  to  inject  self-inflicted  interference 
between  communications  symbols  such  that  a  non-cooperative  receiver  cannot  resolve 
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the  transmitted  communications  symbols.  Two  cooperative  receivers  are  proposed 
that  mitigate  this  self-interference  using  signal  processing  techniques. 

1.1  Problem  Statement 

Because  many  types  of  wireless  communications  are  broadcast  over  a  large  area, 
the  signals  are  highly  vulnerable  to  third-party  interception  and  monitoring.  Protec¬ 
tion  strategies  at  both  the  bit  and  waveform  level  are  desirable  to  prevent  eavesdrop¬ 
ping.  This  can  be  achieved  by  applying  techniques  to  obstruct  non-cooperative  parties 
from  detecting,  intercepting,  and/or  exploiting  the  communications  signal.  The  goal 
of  a  Low  Probability  of  Detection  (LPD)  signal  is  to  prevent  other  parties  from  being 
aware  that  the  signal  is  present.  For  the  purposes  of  this  thesis,  Low  Probability  of 
Interception  (LPI)  refers  to  difficulty  in  identifying  and  extracting  the  symbols  from 
the  waveform.  Finally,  Low  Probability  of  Exploitation  (LPE)  describes  methods  that 
prevent  recovery  of  the  signal’s  information  at  the  bit  level. 

1.2  Research  Focus 

This  thesis  develops  a  modification  to  Differential  Phase  Shift  Keyed  (DPSK) 
modulation  to  achieve  a  waveform  with  LPI  characteristics.  This  is  achieved  by 
inserting  controlled  intersymbol  interference  (ISI)  in  the  transmitted  waveform  that 
inhibits  accurate  symbol  estimation  in  non-cooperative  receivers.  Additional  signal 
processing  in  the  cooperative  receiver  must  be  designed  to  mitigate  the  injected  ISI. 
Symbol  estimation  accuracy  of  the  cooperative  receiver  for  the  LPI  signalling  scheme 
should  be  similar  to  the  optimal  DPSK  signalling  performance.  The  additional  costs 
of  computational  complexity  in  the  cooperative  receiver  and  potential  bandwidth 
expansion  of  the  waveform  are  considered. 

1 . 3  Implications 

The  LPI  signalling  strategy  proposed  in  this  thesis  potentially  provides  an  ad¬ 
ditional  layer  of  information  security  for  wireless  communications  that  can  be  used  in 
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concert  with  existing  LPD/LPI/LPE  protection  schemes.  The  concept  presented  in 
this  research  is  applicable  to  common  modulation  techniques  such  as  Pulse  Amplitude 
Modulation,  Quadrature  Amplitude  Modulation,  and  Phase  Shift  Keyed  Modulation. 

1.4  Preview 

Chapter  [H]  contains  background  information,  outlining  existing  physical  layer 
strategies  for  LPI  communications,  develops  a  model  framework  for  multipath  channel 
effects  and  equalization,  and  surveys  blind  equalization  techniques.  Proposed  signal 
processing  techniques  for  generating  low-probability  of  intercept  waveforms  are  de¬ 
veloped  in  Chapter  III[  Chapter  IV  develops  a  numerical  simulation  of  these  signal 
processing  techniques,  simulation  results  are  presented  in  Chapter  |2J  and  conclusions 
are  summarized  in  Chapter  PVTl 
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II.  Background  and  Literature  Review 


This  chapter  begins  by  reviewing  several  common  existing  strategies  for  creating 
secure  wireless  communications  signals.  LPE  signalling  via  encryption,  LPD 
through  spread  spectrum,  and  LPI  via  frequency  hopping  are  presented.  The  remain¬ 
der  of  the  chapter  discusses  ISI  in  the  context  of  wireless  communications.  The  LPI 
signalling  scheme  developed  in  this  thesis  uses  ISI  to  degrade  non-cooperative  receiver 
performance.  Significant  research  exists  concerning  the  mitigation  of  ISI  caused  by 
multipath  channels.  These  techniques  are  applicable  to  the  signal  processing  design 
of  a  cooperative  receiver  for  the  LPI  scheme  developed  in  Chapter  Hill  Blind  equaliza¬ 
tion  methods,  which  mitigate  ISI  due  to  multipath  with  minimal  a-priori  knowledge 
of  the  signal,  are  also  reviewed  for  potential  application  to  non-cooperative  receivers 
attempting  to  demodulate  the  proposed  LPI  technique. 

2.1  Existing  Secure  Communications  Strategies 

Encryption  is  perhaps  the  most  widely  implemented  form  of  digital  communica¬ 
tions  protection  in  both  the  wireless  and  wired  environments.  Encryption  techniques 
are  widely  studied,  and  a  survey  of  encryption  is  beyond  the  scope  of  this  document. 
It  is  sufficient  to  state  that  the  goal  of  encryption  is  to  prevent  a  third  party  from 
extracting  information  from  an  intercepted  bit  stream.  In  other  words,  encryption 
protects  the  information  in  the  event  that  the  transmitted  bits  are  intercepted.  A 
cooperative  receiver  relics  on  a-priori  knowledge  of  the  decryption  key  to  extract 
the  message  from  the  received  bit  stream.  The  security  of  encryption  relics  on  the 
inability  of  a  third  party  to  determine  the  key,  even  if  the  encryption  technique  is 
known  [21:463] 

Spread  spectrum,  on  the  other  hand,  makes  it  difficult  for  a  third  party  to 
recognize  that  a  communications  signal  is  being  transmitted.  The  modulated  signal 
is  spread  over  a  large  bandwidth  by  mixing  with  a  high  rate  pseudorandom  sequence. 
By  dispersing  the  signal  power  over  a  large  bandwidth,  a  non-cooperative  receiver  is 
less  likely  to  discern  the  presence  of  a  transmitted  signal  from  the  channel  noise  [21:95]. 
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The  pseudorandom  code  used  to  create  a  spread  spectrum  waveform  must  be  known 
by  the  receiver  to  recover  (despread)  the  waveform  to  achieve  to  original  modulated 
communications  signal.  In  many  cases,  a  non-cooperative  receiver  will  not  have  access 
to  this  code,  making  demodulation  difficult  even  if  the  signal  of  interest  is  known  to 
be  present  [15]. 

Frequency  hopping  involves  adjusting  the  carrier  frequency  of  the  modulated 
signal  according  to  a  pseudorandom  sequence.  A  cooperative  receiver  uses  knowledge 
of  the  pseudorandom  sequence  of  carrier  frequencies  to  locate  and  demodulate  the 
signal.  Presumably,  a  non-cooperative  receiver  lacks  access  to  the  hopping  sequence 
and  must  detect  each  successive  carrier  frequency  and  match  the  frequency  in  order 
to  demodulate  the  signal.  To  prevent  signal  interception,  the  hop  rate  must  be  suf¬ 
ficiently  fast  enough  to  prevent  the  non-cooperative  receiver  from  adjusting  to  and 
extracting  information  at  any  carrier  frequency  [21:199]. 

2.2  Intersymbol  Interference  and  Equalization 

Direct  sequence  digital  communications  signals  transmit  information  as  a  series 
of  consecutive  symbols  representing  binary  data.  In  a  basic  signalling  scheme,  symbols 
are  transmitted  in  succession,  with  the  end  of  the  current  symbol  coinciding  exactly 
with  the  beginning  of  the  next  symbol.  The  bandlimited  nature  of  physical  channels, 
however,  limits  the  speed  at  which  the  communications  waveform  can  be  changed. 
This  temporal  spreading  of  the  signal  blurs  the  boundaries  between  communications 
symbols,  creating  ISI  [4,16].  In  wireless  communications,  ISI  is  exacerbated  by  the 
multiple  paths  a  transmitted  signal  may  take  to  reach  a  receiver,  with  each  path 
inducing  different  amplitude  and  phase  (time  delay)  distortions  on  the  signal.  At  any 
single  instant,  the  composite  signal  measured  at  the  receiver  is  the  result  of  multiple 
symbols  that  have  travelled  over  multiple  routes  [16:771]. 

The  composite  effect  of  these  distortions  can  be  modeled  as  a  linear  time- 
invariant  system  characterized  by  the  overall  channel  impulse  response  h(t)  which 
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Figure  2.1:  Baseband  communications  model 


maps  the  transmitted  signal  x(t)  to  the  received  signal  y{t)  according  to  f]2.1jl  [6:70]. 


y(t) 


x  (A)  h(t  —  A)  d\ 


(2.1) 


This  chapter  considers  the  baseband  communications  model  shown  in  Figure  2.1. 
The  transmitted  symbol  sequence  is  represented  by  a  series  of  discrete,  complex  sym¬ 
bols  denoted  x(n).  The  digital  to  analog  converter  (DAC)  block  converts  the  discrete 
series  to  a  continuous  time  ( ct )  impulse  train  with  symbol  impulses  occurring  ev¬ 
ery  Tsym  seconds.  The  transmitted  signal  is  subject  to  the  effective  channel  impulse 
response  h^(t),  which  may  introduce  ISI,  and  additive  noise  b^ct\t).  The  effective 
channel  impulse  response  represents  pulse  shaping  in  the  transmitter  /lrf)(t),  the  phys¬ 
ical  channel  impulse  response  h^ys(t),  and  a  matched  filter  in  the  receiver  t). 

The  superscript  ( ct )  denotes  continuous  time  signals  in  this  thesis.  The  received  sig¬ 
nal  after  matched  filtering  is  y ^  (t)  and  is  sampled  an  integer  number  of  times  O  per 
symbol  period  Tsym. 

The  goal  of  equalization  is  to  mitigate  or  even  remove  the  ISI  via  signal  process¬ 
ing  within  the  communications  receiver.  Filtering,  a  common  approach  to  equaliza¬ 
tion,  attempts  to  modify  the  received  signal  in  a  manner  that  isolates  the  individual 
symbols,  essentially  inverting  the  effect  of  the  channel  [16:150].  The  digital  transverse 


filter  is  a  linear  time-invariant  system  defined  by  its  discrete  impulse  response  w' 
which  acts  on  the  received  signal  samples  y'(k )  to  produce  an  estimate  of  the  trans¬ 
mitted  signal  (x(n)  ~  x(n  —  5))  as  shown  in  Figure  12.11  The  “prime”  designation 
denotes  oversampled  signals  with  Tsym/0  seconds  between  samples.  The  transverse 
filtering  operation  is  described  mathematically  by  the  discrete  convolution  sum: 

M—l 

x\k)  =  w\k)  *  y\k)  =  ^  w'(m)y'(k  —  m)  .  (2.2) 

m= 0 

The  above  model  assumes  the  channel  response  is  static.  In  many  physical 
systems,  however,  the  channel  impulse  response  varies  over  time  [16:946].  For  static 
channels,  preset  equalization  in  which  the  filter  taps  are  constant  is  sufficient.  For 
slowly  time- varying  channels,  an  adaptive  equalizer  is  required.  An  adaptive  equalizer 
may  update  the  filter  taps  continuously  or  periodically.  In  either  case,  an  adaptive 
filter  is  characterized  by  a  convergence  rate.  When  the  channel  variation  rate  exceeds 
the  convergence  rate  of  the  filter,  symbol  recovery  is  degraded  [16:158].  The  LPI 
signal  design  in  Chapter  [III]  attempts  to  vary  the  effective  channel  impulse  response 
beyond  the  filter  convergence  rate  to  prevent  eavesdropping. 

Of  particular  interest  to  non-cooperative  receiver  design  are  blind,  adaptive 
equalization  methods.  These  methods  are  designed  to  undo  the  effects  of  the  linear, 
time- invariant  channel  using  only  limited  knowledge  of  the  channel  and  the  statistics 
of  the  source  signal  [4:684],  A  successful  blind  equalizer  allows  a  non-cooperative 
receiver  to  minimize  ISI  and  extract  communications  signals  with  little  or  no  prior 
knowledge  of  the  information  content. 

2.3  Single/ Multichannel  Communications  Model  Equivalence 

When  the  sampling  period  is  an  integer  multiple  of  the  symbol  period  in  the 
single  channel  model  (Figure  I2.ll).  the  system  can  be  equivalently  represented  as  a 
multiple  channel  system  with  discrete  signals  and  impulse  responses  at  a  sampling 
frequency  of  1/Tsym  as  shown  in  Figure  [2i21  This  section  derives  the  eqiiivalent  multi- 
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Noise 
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channel  model  and  presents  a  block  matrix  notation  for  the  system.  The  equivalence 
derivation  is  based  on  [5]. 

2.3.1  Model  Equivalence.  In  this  thesis,  the  effective  channel  impulse  re¬ 
sponse  is  a  combination  of  the  transmitter  pulse  shape,  the  physical  channel  model, 
and  the  receiver  matched  filter.  The  continuous  time  effective  channel  impulse  re¬ 
sponse  h^jj  (t)  can  be  expressed  as 

»  /'“>*(-()  ,  (2.3) 


where 


/(rf)(t)  =  Pulse  shape  (and  matched  filter)  , 
h physi t)  =  Physical  channel  impulse  response. 


All  pulse  shapes  considered  in  this  thesis  are  real  and  symmetric,  so  the  pulse  shape 
and  matched  filter  are  identical  £)).  Note  that  continuous  time 

signals  are  explicitly  denoted  by  the  superscript  (ct). 
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In  the  baseband  communications  model,  the  N  transmitted  symbols  are  repre¬ 
sented  by  a  series  of  complex,  discrete  terms  x(ri).  The  continuous  time  signal  after 
matched  filtering  in  the  receiver  can  be  written  as 

N- 1 

y^it)  =  ^2  x(n)h%pf  {t  -  nTsym  -  tQ)  +  6(ct)(t)  .  (2.4) 

n= 0 


The  additive  noise  signal  b(ct>  (t)  is  the  channel  noise  passed  through  the  matched 
filter.  Unless  otherwise  stated,  the  channel  noise  is  assumed  to  be  additive  white 
Gaussian  noise  It  follows  that  the  resulting  noise  signal  is  described  by 


(2.5) 


The  oversampling  factor  O  defines  the  number  of  samples  per  symbol  period 
Tsym,  resulting  in  a  sampling  rate  of  0/Tsym  samples  per  second.  In  this  way,  the 
sampled  matched  filter  output  can  be  written  as 


The  synchronization  of  the  sampling  is  determined  by  t0.  The  sampled  signals 
in  (2.6)  have  obvious  discrete  analogs  defined  in  (277).  The  “prime”  designation  is  used 
to  indicate  an  oversampled  discrete  series  with  Tsym/0  seconds  between  samples. 

V'(k)  = !,'«) 

h'(k)  =  (T=  -  t„) 

b'(k)  =  &<“*>  {kT‘J"  '\  (2.7) 
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After  sampling  the  matched  filter  output  in  the  baseband  model,  the  equaliza¬ 
tion  filter  with  M  taps  spaced  Tsym/0  seconds  apart  is  applied.  The  resulting  signal 
x'{k)  is  given  by 

M— 1 

x'(k)  =  ^  w'(m)y ^ 

m= 0 

The  approximated  symbols  are  extracted  by  downsampling  by  decimating  by  O: 
x(n)  =  x'  ( 0(n  +  1)  —  1)  , 

M- 1 

=  ^  w\m)y ^ 

m= 0 

=  x:  w'(m)y(a)  (nTnm  +  {— •  (2.9) 

m= 0  '  ' 


( 0(n  +  1) 


m) 


■  sym 

IT 


ik-rn)^ 


(2.8) 


The  multichannel  model  can  be  extracted  from  (j2.9j)  by  restructuring  as  O 
summations  of  decimated  series  w'(m )  and  y'(k ): 

(M/O)  — l  /  \ 

x[n)  =  w\Om)y{ct)  Un-m)  Tsym  +  - -  sym  J 

m= 0  '  ' 

(M/O)- 1  /  (O  —  9  IT  \ 

+  w'{Om  +  l)//(rf)  (n-m)  J 

m=0  '  ' 

(M/O)- 1 

+  ^  w\Om  +  O  -  l)y(ct)  ((n  -  m)  Tsym)  , 

m=0 

(M/O)  — 1 

=  ^  w'(Om)?/  ((n  -  m)0  +  O  -  1) 

ra=0 
(M/O)— 1 

+  ^  w'(Oro  +  1)?/  ((n  -  m)0  +  O  -  2) 

m= 0 

(M/O)  — 1 

+  ^  +  0-1)?/  ((n  -  m)0)  .  (2.10) 

m=0 
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Defining  sample  spaced  subchannel  variables  in  terms  of  the  oversampled  single 
channel  variables  (12.1211  allows  the  estimated  symbols  to  be  written  as  O  discrete 
convolutions: 

x(n)  —  (w^\n)  *  y(°\n))  +  (w^2\n)  *  + . . .  +  (w(°\n)  *  y^^n))  .  (2.11) 


b{i){n)  =  b'  (nO  +  i-  1) 
y^l\n)  =  y'  ( nO  +  i  —  1) 

h®{n)  =  h(nO  +  i- 1)  (2.12) 


Finally,  the  subchannel  siganls  y^\n)  are  derived  in  terms  of  the  subchannel 
impulse  responses: 


JV-l 


y^\m)  =  x(n)h 


M  ((°m  +  i  1)7^=  -  „ T,vm  -tc  )  +  b «>(m)  , 


71— 0 

N—l 


eff 


^  x(n)h^j  (  (m  -  n)Tsym  + 


O 


71—0 


(*  -  1)T« 


sym 


o 


-  to  )  +  few(m) 


=  x(m)  *  hSl\m)  +  b^\m)  . 


(2.13) 


So  hb),  b^(n),  y®(n),  and  are  now  defined  in  terms  of  their  single  channel 
analogs.  With  these  definitions,  (l2.ll)  and  (2.1311  map  the  system  input  x(n)  to  the 
output  x[n)  and  the  multichannel  model  is  described  completely. 

2.3.2  Convolution  Matrix.  The  discrete  linear  convolution  operation  can  be 
expressed  in  matrix  form.  In  general,  recasting  the  convolution  as  a  matrix  opera¬ 
tion  allows  for  faster  convolution  computations.  In  the  context  of  communications 
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channel  modeling  and  equalization  filtering,  properties  of  the  convolution  matrix  cor¬ 
responding  to  the  effective  channel  impulse  response  provide  insight  into  equalization 
performance  (see  Section  12.4.3).  The  following  block  notation  is  derived  from  [5] 
and  [11]. 

To  begin,  define  the  subchannel  impulse  response  h^\n)  and  equalization  filters 
w^\n)  sequences  as  vectors.  Each  impulse  response  vector  has  length  L  and  each 
equalization  filter  has  length  M/O. 

h(<)  =  [/i«(0),  hw(  1),  ...  h{i){L  -  1)]T  (2.14) 

w«  =  [ww(0),  ww(l),  •••  w{i){M/0  -  1)]T  (2.15) 


The  subchannel  equalization  hlters  can  be  stacked  columnwise  to  make  a  single  vector 
w (n)  according  to  (2.161).  which  should  not  be  confused  with  the  oversampled  filter 
tap  vector  w'  in  Figure  j2.lt 


w  (n) 


(2.16) 


The  subchannel  matched  filter  output  y^\n)  can  be  expressed  as  a  vector  by 
recasting  all  of  the  discrete  sequences  in  terms  of  vectors  or  convolution  matrices.  To 
proceed,  let  y  w(n)  be  a  vector  of  N  matched  filter  samples  such  that 

yW(n)  =  [y^\n),  y{l\n  -  1),  yM(n  -  2),  ...  y[l)(n  -  N  +  1)]T  .  (2.17) 

Similarly,  x(n)  and  b^(n)  are  defined  according  to 

x(n)  =  [x(n),  x(n  —  1),  x(n  —  2),  ...  x(n  —  N  —  L  +  2)]T  ,  (2.18) 

b«(n)  =  [b^(n),  b{i\n  -  1),  b{i\n  -  2),  ...  6(i)(n  —  +  1)]T  .  (2.19) 
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The  subchannel  convolution  matrix  for  the  ith  subchannel  acts  on  (N  + 
L—  1)  transmitted  symbols  in  x^(n)  to  give  N  samples  in  y^(n)  (|2.2Q[) .  The  length 
of  the  effective  channel  impulse  response  as  an  integer  number  of  symbol  periods  is  L. 
is  a  block  diagonal  matrix  with  on  the  diagonal  (2.21). 


yw(n) 


H^x(n)  +  b^(n) 


(2.20) 


N+L 


h®(  0)  ••• 

h^(L-l)  0  . 

0 

rO)  _ 

nN  ~ 

0  /i«(  0) 

1 

O 

0 

> 

0 

0  h®(0)  •• 

•  h®(L-  1) 

> 

>N 


(2.21) 


Stacking  all  of  the  subchannels  columnwise,  the  entire  multichannel  matched 
filter  output  y  (n)  can  be  expressed  in  terms  of  the  block  convolution  matrix  containing 
the  effective  channel  impulse  response  of  each  subchannel: 


y  (n)  = 


y  (1\n) 


y(o)(n) 


H 


(i) 

N 


H 


(O) 

N 


x  (n 


bW(n) 

b^(n) 


=  Hjvx(n)  +  b(n) 


(2.22) 


Calculating  the  estimated  symbol  vector  x(n)  of  length  N  requires  a  matched 
filter  output  vector  y (n)  of  length  0(N  +  M/O  —  1).  Each  subchannel  equalization 
filter  can  be  cast  as  a  convolution  matrix  wjy  with  dimensionality  ( M/O  + 
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N)  x  N.  Concatenating  these  convolution  matrices  row-wise  (12.231)  produces  the 
estimated  symbol  vector  solution  in  (12.241). 


wN= 


(2.23) 

(2.24) 


x(n)  =  Wjvy(n) 


2.4  Equalization  Techniques 

To  better  understand  the  major  equalization  algorithm  families,  a  classification 
structure  is  desired.  The  equalization  algorithms  discussed  in  this  section  are  shown 
in  Figure  12731  organized  according  to  the  hierarchy  presented  in  Haykin’s  Adaptive  Fil¬ 
ter  Theory  [4] .  This  survey  contrasts  several  key  algorithm  characteristics  relevant  to 
non-cooperative  receiver  operation  as  discussed  in  Section  12.4.11  After  defining  these 
algorithm  traits,  the  algorithms  of  interest  are  summarized.  Because  the  behavior  of 
blind  equalization  algorithms  is  often  described  relative  to  least  mean  square  (LMS) 
adaptive  filtering  based  on  a  known  training  signal,  an  overview  of  LMS  is  presented 
first  [4:231].  The  constant  modulus  and  Sato  algorithms,  which  share  a  similar  struc¬ 
ture  with  the  LMS  algorithm,  are  presented  next.  A  block  subspace  decomposition 
algorithm  based  on  eigenvalue  decomposition  of  the  received  signal’s  autocorrelation 
matrix  is  described  in  detail.  Finally,  spectrum  fitting  and  Shalvi  and  Weinstein’s 
higher-order  cumulants  approaches  are  reviewed.  Special  attention  is  given  to  the 
statistical  assumptions  made  by  each  algorithm  and  the  performance  under  non-ideal 
conditions  for  possible  LP1  signalling  techniques. 

2-4-1  Equalization  Algorithm  Key  Characteristics.  Haykin  groups  blind 
algorithms  into  two  broad  categories:  algorithms  based  on  the  estimation  of  a  signal’s 
second  order  statistics  (SOS)  vice  algorithms  requiring  higher-order  statistics  (HOS) 
to  be  estimated.  Because  second  order  statistics  can  be  estimated  accurately  with 
far  fewer  observations  than  higher-order  statistics,  SOS  based  algorithms  typically 
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Figure  2.3:  Hierarchy  of  equalization  algorithms 


2-12 


converge  significantly  faster  than  those  leveraging  HOS  [4],  Convergence  speed  can 
be  quite  important  in  a  wireless  environment  with  a  dynamic  channel  response.  The 
ability  of  the  receiver  to  extract  meaningful  information  often  hinges  on  whether 
the  equalizer  can  track  the  changing  impulse  response  of  the  channel  and  invert  the 
multipath/fading  effect  [24], 

Equalizer  sampling  rate  is  another  characteristic  useful  for  classification.  Gen¬ 
erally,  an  algorithm  either  uses  a  single  sample  from  each  communications  symbol, 
or  multiple  samples  (note  that  this  assumes  the  symbol  rate  is  known  by  the  re¬ 
ceiver).  Algorithms  utilizing  multiple  samples  per  transmitted  symbol  are  considered 
fractionally  spaced  equalizers  (FSE),  and  generally  perform  better  than  equivalent 
sample  spaced  or  baud  spaced  equalization  (BSE)  methods  [16:161].  Fractional  sam¬ 
pling  techniques  often  converge  faster  and  with  greater  accuracy  than  their  sample 
spaced  analogs  [4:727]. 

The  ability  of  a  blind  equalization  algorithm  to  adapt  in  real  time  is  useful 
to  consider,  particularly  when  considering  applications  in  wireless  systems  which  are 
often  subject  to  time-varying  channel  behavior.  Algorithms  that  update  the  equal¬ 
ization  filter  on  a  symbol-by-symbol  basis  are  considered  to  be  continuously  adaptive. 
Algorithms  which  process  a  block  of  observed  symbols  transmitted  under  a  single 
channel  condition  to  identify  an  equalization  filter  are  labeled  periodically  adaptive 
for  the  purposes  of  this  thesis. 

An  equalization  method  is  considered  direct  if  the  output  of  the  algorithm  is 
the  tap  weights  of  the  linear  equalization  filter  w.  For  example,  the  equalization 
filter  is  estimated  many  times  in  some  iterative  algorithms.  The  indirect  equaliza¬ 
tion  approach  first  estimates  the  channel  impulse  response  h  before  solving  for  the 
equalization  filter  taps. 

2.4.2  Least  Mean  Square  (LMS)  Algorithm.  The  LMS  algorithm  applied  to 
equalization  is  an  adaptive  method  of  updating  the  tap  weights  of  a  transverse  filter 
based  on  the  difference  between  the  filter  output  and  the  known  source  signal  originally 
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Figure  2.4:  Block  diagram  of  the  least  mean  square  algorithm 

transmitted.  The  process  is  a  stochastic  extension  of  deterministic  steepest  descent 
algorithms  that  use  the  gradient  of  an  error  term  to  determine  how  the  tap  weights 
should  be  updated  to  decrease  the  error  term.  Instead  of  explicitly  computing  the 
gradient,  the  LMS  algorithm  estimates  the  gradient  based  on  the  known  transmitted 
signal,  received  signal,  and  current  filter  tap  weights.  Figure  12.41  depicts  a  receiver 
implementation  of  the  LMS  algorithm  [4:231].  The  equalization  filter  is  updated  after 
processing  each  received  symbol,  giving  a  new  set  of  tap  coefficients  for  w  at  each 
symbol  index  n. 

Both  the  method  of  steepest  descent  and  LMS  algorithms  use  a  cost  function 
defined  by  the  expectation  of  the  square  of  the  error  term  e(n).  The  resulting  cost 
function  can  be  written  in  terms  of  the  statistics  of  the  transmitted  and  received 
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signals  [18:287].  Using  the  multichannel  model  notation,  the  cost  function  is 


J(n)  =  E  [e*(n)e(n)\  , 

=  E  [(x(n  —  5)  —  wr(n)y(n))*  (x(n  —  5)  —  wT(n)y(n))]  , 

=  E  [x*(n  —  S)x(n  —  5)]  —  2Re  <  E  \x(n  —  S)yH(n)]  w*(n)  > 

\  Ply  ) 

+  w T(n)  E  [y{n)yH{n)]  w*(n)  , 

" - . - ' 

R  y 

=  E  [x*(n  —  5)x{n  —  <5)]  —  2Re  {p^w*(n)}  +  wT  (n)Ryw*(n)  .  (2.25) 


Taking  the  gradient  of  the  cost  function  with  respect  to  the  filter  w  gives 

Vw  J(n)  =  -2p^y  +  2wr(u)Rj/  .  (2.26) 

Note  that  the  expectations  have  been  rewritten  as  the  correlation  matrix  R;/  of  the 
received  signal  vector  y (n)  and  the  cross  correlation  vector  pxy  of  the  received  signal 
y(n)  with  the  current  symbol  from  the  known  training  signal  [18:287].  Convergence 
of  the  steepest  descent  method  depends  on  the  autocorrelation  matrix  and  cross¬ 
correlation  vector  remaining  constant  at  all  time  steps  n,  thus  the  transmitted  symbol 
sequence  x{n)  is  assumed  to  be  wide-sense  stationary  (WSS)  and  the  received  signal 
after  sampling  is  either  WSS  if  y'[k)  is  sample  spaced  or  cyclostationary  if  y'[k)  is 
fractionally  spaced.  In  practice,  LMS  can  be  effective  for  time-varying  channels  if  the 
convergence  speed  is  fast  relative  to  the  changing  channel  response,  even  though  this 
violates  the  stationary  assumption  for  y{n)  [4:203-205]. 

In  LMS,  the  gradient  of  the  cost  function  is  no  longer  a  deterministic  function  of 
filter  coefficients.  Instead,  the  gradient  at  each  time  step  is  estimated  stochastically 
using  a  different  estimate  of  the  autocorrelation  and  cross-correlation  matrices  at  each 
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time  step  according  to 


VwJ(n)  =  — 2p T{n)  +  2w(n)TRy(n)  , 


(2.27) 


where 


&v(n)  =  y(n)y H(ri)  , 

P xy(n)  =  y *(n)x(n  -  S )  . 

Note  that  the  length  of  the  received  vector  must  match  the  length  of  the  filter.  The 
filter  length,  in  turn,  needs  to  be  large  enough  to  approximate  the  inverse  channel 
response.  Using  three  times  the  length  of  the  finite  channel  impulse  response  model 
to  determine  the  filter  length  is  typical  [4:236]. 

The  estimated  gradient  of  the  cost  function  produces  the  filter  update  function 
given  by 

w(n  +  1)  =  w (n)  —  i-i  ^VJ(n) j  , 

=  w(n)  +  fiy*(n)  (x(n)  —  yT(n)w(n))  , 

=  w(n)  +  /iy*(n)e(n)  .  (2.28) 

As  implied  by  the  function,  the  adaptive  controller  in  the  LMS  algorithm  updates 
the  filter  after  each  received  symbol  [4:236].  The  parameter  /i  is  the  step  size  that 
controls  how  much  the  filter  weights  are  updated  in  the  direction  indicated  by  VJ(n). 
The  step  size  must  be  small  enough  to  prevent  the  filter  coefficients  from  jumping 
too  far  past  the  error  minima  with  each  iteration  and  diverging.  For  values  of  /z  that 
converge,  however,  smaller  step  sizes  result  in  slower  convergence.  Although  the  step 
size  is  somewhat  arbitrary,  the  general  convergence  criteria  calls  for  /i  to  be  small 
relative  to  the  inverse  of  the  largest  eigenvalue  of  Ry  [4:288]. 
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Because  the  symbol-by-symbol  update  of  the  filter  coefficients  is  based  on  an 
estimate  of  the  gradient  that  varies  with  every  iteration,  the  filter  never  converges  to 
the  exact  equalization  solution  that  inverts  the  channel  effects  completely.  Instead,  the 
estimate  will  continue  to  fluctuate  around  the  optimal  solution  after  convergence.  The 
average  residual  mean  square  error,  or  misadjustment,  after  convergence  is  directly 
proportional  to  the  step  size.  Thus,  step  size  selection  must  balance  convergence 
speed  and  steady  state  error.  This  has  driven  investigation  of  LMS  algorithms  with 
variable  step  sizes  [7]. 


2-4-3  Constant  Modulus  Algorithm  (CMA).  The  constant  modulus  algo¬ 
rithm  (CMA)  is  one  of  the  most  common  blind  equalization  techniques  [5].  The 
algorithm  is  functionally  similar  to  the  LMS  algorithm,  sharing  the  same  adaptive 
update  method  and  modest  processing  requirements  for  implementation.  Both  sym¬ 
bol  spaced  and  fractionally  spaced  versions  have  been  studied  extensively.  CMA  is 
a  part  of  a  larger  class  of  Bussgang  blind  adaptive  algorithms  that  directly  equalize 
the  channel  based  on  higher-order  statistics  [4].  The  fractionally  spaced  version  has 
the  advantage  of  converging  to  the  same  equalization  solution  as  the  trained  LMS 
algorithm  (perfect  equalization)  under  noiseless  conditions  for  most  channels  [5].  The 
criteria  for  convergence  to  the  LMS  solution  are  presented  below.  Sample  spaced 
CMA,  on  the  other  hand,  has  been  shown  to  require  impractical  constraints  such  as 
infinite  filter  length  in  order  to  achieve  perfect  equalization  [2], 

Intuitively,  CMA  works  by  minimizing  an  error  term  based  on  the  difference 
between  the  expected  symbol  amplitude  and  the  sampled  magnitude  of  the  symbol: 


e{n) 


E[|s(n)|4] 
E  [\x(n)\2] 


x(ri)\ 2  . 


(2.29) 


The  error  term  can  be  minimized  using  the  same  update  strategy  as  the  LMS  error 
minimization  technique  described  in  Section  12.4.21  Error  minimization  operates  on 
the  square  of  the  error  term,  which  itself  contains  the  symbol  estimate  squared  \x(n)\2, 
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so  convergence  relies  implicitly  on  the  estimation  of  fourth  order  statistics  via  time 
averaging.  This  HOS  dependence  dictates  the  low  convergence  speed  of  the  CMA 
algorithm  relative  to  SOS  based  methods  [4:730]. 

Clearly,  if  the  source  signal  is  of  constant  modulus,  the  error  term  approaches 
zero  as  the  filter  approaches  convergence.  Although  not  obvious,  CMA  converges 
even  when  the  source  signal  does  not  have  a  constant  modulus.  In  this  case,  however, 
the  minima  of  the  cost  surface  no  longer  correspond  to  zero  cost.  The  non-zero 
steady  state  cost  of  the  minima  have  the  effect  of  “flattening”  the  shape  of  the  cost 
surface,  resulting  in  slower  convergence.  This  “flattening”  phenomena  becomes  more 
pronounced  as  the  kurtosis  of  the  source  signal  increases  [5:1935]. 

As  with  the  LMS  algorithm,  the  step  size  /i  must  strike  a  balance  between 
convergence  speed  and  steady  state  mean  square  error,  making  variable  step  size 
strategies  desirable  [23].  For  example,  one  such  method  performs  a  second  cost  mini¬ 
mization  with  respect  to  the  step  size  at  each  iteration,  searching  the  cost  surface  in 
the  direction  of  steepest  descent  as  calculated  by  the  gradient  estimation.  This  mini¬ 
mization  determines  the  “optimal”  step  size  for  the  current  iteration.  This  approach 
has  the  additional  advantage  of  reducing  the  likelihood  of  converging  to  local  minima 
for  multimodal  cost  surfaces  [26]. 

Like  LMS,  CMA  assumes  a  wide-sense  stationary  or  cyclostationary  received 
signal  (depending  on  the  sampling  rate  of  the  receiver).  Under  time- varying  channel 
conditions,  the  received  signal  statistics  will  fluctuate,  violating  the  WSS  assumption. 
However,  given  sufficiently  slow  channel  variation,  the  signal  statistics  may  still  be 
stable  enough  for  the  time- averaging  implicit  in  CMA  to  equalize  the  channel  and 
track  the  variation  in  cost  surface  minima  [4:688]. 

In  the  case  of  fractionally  spaced  CMA,  the  following  four  conditions  guarantee 
that  the  cost  surface  minima  correspond  to  Liter  coefficients  that  perfectly  equalize  the 
channel  response.  First,  the  channel  must  be  noiseless.  Second,  the  (possibly  complex) 
source  signal  must  be  composed  of  a  series  of  independent,  identical  symbols  from  a 
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circularly  symmetric  constellation.  Third,  the  kurtosis  of  the  symbol  distribution 
must  be  less  than  the  kurtosis  of  a  Gaussian  distribution  [4:727].  Lastly,  the  channel 
convolution  matrix  as  defined  by  Tong  et  al.  must  be  full  column  rank  [20].  The 
combined  channel-filter  impulse  response  for  an  FSE  is  presented  in  terms  of  the 
effective  channel  convolution  matrix  in  Section  12.3.21  Notably,  the  second  and  third 
requirements  for  perfect  equalization  are  properties  of  the  source  signal,  suggesting 
that  these  features  may  be  exploitable  in  the  design  of  LPI  communications  systems, 
although  these  assumptions  are  not  exploited  in  this  thesis. 

The  full  column  rank  requirement  of  the  channel  convolution  matrix  leads  to 
a  minimum  number  of  filter  taps  required  for  perfect  equalization.  For  a  matrix  to 
be  full  column  rank,  it  must  have  at  least  as  many  rows  as  columns.  The  number  of 
required  taps  is  a  function  of  the  length  of  the  channel  response  and  the  number  of 
samples  during  each  symbol  period  according  to  (12.30b  where  M  is  the  total  number 
of  filter  tapes,  either  arranged  in  a  single,  fractionally  spaced  vector  in  the  single 
receiver  model,  or  as  O  symbol  spaced  vectors  of  length  M/O  in  the  multichannel 
model.  L  is  the  number  of  symbol  periods  over  which  the  channel  impulse  response 
is  non-zero,  and  O  is  the  number  of  samples  taken  during  each  symbol  period  in  the 
single  channel  model,  or  the  number  of  antennas  in  the  multichannel  model  [5:1930]. 

m>l  +  M-  i 

M  ~  0{(L)-\)  ~  L  (2.30) 


Another  extension  of  the  full  column  rank  condition  is  the  requirement  that 
the  impulse  response  of  each  channel  (in  the  multichannel  system  model)  be  linearly 
independent  of  other  channel  impulse  responses.  In  other  words,  the  z-transform  of 
the  channel  responses  must  share  no  common  zeros  [4:693]. 
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Many  variations  of  CMA  have  been  developed,  typically  with  the  goal  of  ad¬ 
dressing  specific  CMA  shortcomings.  Vector  CMA  [25]  is  of  particular  interest  in 
relation  to  the  LPI  strategy  proposed  in  Chapter  11111  Vector  CMA  suggests  a  means 
of  equalizing  source  signals  with  high  kurtosis  such  as  QAM  constellations  produced 
by  shell  mapping.  Vector  CMA  defines  an  error  term  based  on  a  vector  of  successive 
symbols  in  an  attempt  to  create  a  statistic  with  sub-Gaussian  kurtosis.  In  the  case 
of  shell  mapping,  the  encoding  structure  is  such  that  a  properly  chosen  symbol  block 
length  produces  uniformly  distributed  error  statistics  that  can  be  minimized  via  the 
cost  surface  gradient  descent  [22,25].  This  process  is  similar  to  the  bank  of  LMS 
equalizers  discussed  as  a  cooperative  receiver  structure  for  the  LPI  signal  design. 

2-4-4  Suto  Algorithm.  Like  the  constant  modulus  algorithm,  Sato’s  algo¬ 
rithm  is  an  adaptive  algorithm  based  on  LMS.  As  in  LMS  and  CMA,  the  Sato  algo¬ 
rithm  directly  equalizes  the  received  signal  by  iteratively  adjusting  the  equalization 
filter  coefficients.  As  originally  proposed  in  [13],  the  Sato  algorithm  is  developed  in 
the  context  of  multilevel  pulse  amplitude  modulation  (PAM).  The  algorithm  simplifies 
equalization  by  dividing  the  PAM  detection  space  into  two  equally  likely  regions,  and 
then  modelling  the  signal  as  binary  with  the  remaining  levels  considered  as  additive 
noise.  In  a  zero-mean  PAM  signal,  this  is  equivalent  to  considering  only  the  sign  of 
the  signal  [13] .  The  resulting  error  term  for  a  received  symbol  is  given  by 

e(n)  =  x(n)  —  [sgn  x(n)\  ■  E  [|x(n)|]  .  (2-31) 

The  Sato  algorithm  is  subject  to  the  same  channel  impulse  response  constraints 
as  CMA.  Local  minima  in  the  cost  surface  are  also  a  difficulty,  particularly  in  noisy 
environments.  With  no  particular  advantages  over  CMA  and  an  infinite  equalizer 
length  requirement  for  perfect  equalization  [4:722],  the  Sato  algorithm  is  not  pursued 
as  a  non-cooperative  receiver  architecture  in  this  thesis. 
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2-4-5  Block  Subspace  Decomposition  Algorithm.  Subspace  decomposition  is 
a  channel  identification  algorithm  (indirect  equalization)  for  oversampled  signals  with 
cyclostationary  statistics.  The  primary  advantage  of  the  subspace  decomposition 
method  over  CMA  is  that  it  relies  on  second-order  statistics,  which  can  be  estimated 
well  with  relatively  few  samples  compared  to  higher-order  statistics.  The  algorithm  is 
not  continually  adaptive  in  the  sense  that  the  channel  is  not  estimate  is  not  iteratively 
refined  on  a  symbol  by  symbol  basis.  Instead,  a  block  of  consecutive  samples  must  be 
collected.  Conceivably,  a  sufficiently  slow  time-varying  channel  could  be  repeatedly 
estimated  using  a  sliding  block  of  samples,  making  it  a  candidate  for  a  periodically 
adaptive  implementation.  The  significant  complexity  of  subspace  decomposition  is  a 
limiting  factor.  Also,  as  a  fractionally  spaced  algorithm,  subspace  decomposition  is 
subject  to  the  same  full  column  rank  requirement  on  the  channel  convolution  matrix 
as  CMA  as  discussed  in  Section  [2.4.31  [4:689-701]. 

The  subspace  decomposition  method  considers  the  correlation  matrices  of  the 
transmitted  signal  and  the  received  signal.  Interestingly,  the  theoretical  development 
of  the  subspace  method  includes  an  additive  noise  source  which  is  restricted  to  be  ban- 
dlimited,  stationary,  and  independent  from  the  transmitted  signal  without  disturbing 
the  ability  of  the  algorithm  to  approach  perfect  equalization  [11,20].  Specifically,  the 
correlation  matrix  of  the  received  signal  in  terms  of  the  correlations  of  the  transmitted 
and  noise  vectors  is  derived  in  (2.32).  where  H  is  the  channel  convolution  matrix  as 
derived  in  Section  [2.3.2  [4:694],  Note  that  the  final  line  assumes  that  the  noise  vector 
is  composed  of  independent,  identically  distributed  random  variables  with  variance 


R,  =  E  [yyH] 


=  E 


\H 


(Hx  +  u)  (Hx  +  v 
=  E  [Hxx^Hff]  +  E  [uuH] 
=  HRxHH  +  R„ 

=  HR(H^  +  cr2I 


(2.32) 
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Rewriting  the  correlation  matrix  of  the  received  signal  in  terms  of  eigenvalues 
and  eigenvectors  using  Mercer’s  theorem  (also  called  the  spectral  theorem)  produces 

OJV-l  ON- 1 

'y  ^  =  y  A^q^q^j.  +  <7  I  ,  (2.33) 

k= 0  k= 0 

where  q^*,  and  Xy ^  are  the  eigenvectors  and  eigenvalues  of  the  composite  received 
signal  correlation  matrix  Ry,  q,  ^  and  A Z)k  are  the  eigenvectors  and  eigenvalues  of  the 
noiseless  transmitted  signal’s  correlation  matrix  after  channel  effects,  O  is  the  number 
of  received  samples  per  transmitted  symbol,  and  N  is  the  total  number  of  received 
symbols  [4:813].  Because  each  correlation  matrix  R  is  positive  semi-definite,  the 
eigenvalues  of  each  term  in  (2.33)  are  all  nonnegative.  Although  each  term  produces 
a  square  correlation  matrix  of  size  ON,  the  correlation  matrix  of  the  transmitted 
signal  HRj.Hfl  has  a  rank  of  only  (N  +  L ),  where  L  is  the  channel  length  in  terms  of 
symbol  periods  (effective  channel  order).  Thus.  HRxHfl  has  at  most  (N  +  L)  positive 
eigenvalues  and  ON  —  (N  +  L)  eigenvalues  equal  to  zero.  This  naturally  divides  vector 
space  of  Ry  into  two  subspaces:  the  eigenvectors  corresponding  to  the  (N  +  L)  largest 
eigenvalues  span  the  signal  subspace  and  the  eigenvectors  corresponding  to  the  rest 
of  the  eigenvalues  (each  with  magnitude  cr2)  span  the  noise  subspace  [11]. 

Defining  the  eigenvectors  spanning  the  noise  subspace  as  leads  to  (2.34). 

Substituting  Ry  as  defined  in  (2.32)  into  (2.34)  leads  to  (2.35)  [4:695].  This  result  can 

also  be  reached  directly  by  stating  that  the  noise  subspace  spanned  by  g*,  is  orthogonal 
to  the  subspace  defined  by  the  columns  of  the  convolution  matrix  H.  Moulines,  et 
al.  showed  that  solving  (2.35)  for  H  with  the  restriction  H  ^  0  produces  a  unique 
solution  (with  an  arbitrary  scale  factor)  [11], 

Rygfc  =  cr2gfc,  k  =  N  +  L  +  1, . . . ,  ON  (2.34) 

H^gfc  =  0,  k  =  N  +  L  +  l,...,ON  (2.35) 
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To  cleanly  extract  the  channel  coefficients  h,  (J2.35j)  is  often  restructured.  In  (12.361) , 
each  eigenvector  has  been  recast  into  a  convolution  matrix  G k  following  the  process 
in  Section  2.3.21 

Gf  h  =  0  (2.36) 


Furthermore,  (j2.36j)  can  be  expressed  as  a  scalar  by  setting  the  square  norm  of  G^h 
to  zero: 

||Gfh||2  =  hJ7G^.Gfh  =  0  .  (2.37) 

Finally,  a  cost  function  can  be  defined  in  terms  of  (2.37).  In  this  cost  function  (2.38). 
the  eigenvector  convolution  matrix  is  replaced  by  an  estimate  computed  from  the  block 
of  received  samples.  The  cost  function  can  then  be  minimized  under  the  constraint 
that  h  ^  0  [4:697], 

(ON—N—L—l  \ 

£  GtG?  h  (2.38) 

k= 0  ) 

After  identifying  the  channel  estimate,  the  estimated  channel  convolution  matrix 
H  can  be  constructed  directly.  An  estimate  of  the  transmitted  signal  can  be  obtained 
by  taking  the  pseudoinverse  of  H: 


x  =  HV=  (h*h)  . 


(2.39) 


Alternatively,  using  the  singular  value  decomposition  of  Ry  to  extract  the  eigenvalues 
leads  to  a  transmitted  signal  estimation  method  that  does  not  require  the  additional 
pseudoinverse  computation  [20]. 

In  the  channel  identification  development  above,  several  assumptions  have  been 
made.  As  mentioned  previously,  the  received  signal  must  be  cyclostationary  in  order  to 
estimate  the  correlation  matrix  from  via  time-averaging  [20].  In  addition,  the  channel 
noise  was  assumed  to  be  white  (noise  samples  are  uncorrelated).  The  variance  of 
the  noise  was  assumed  to  be  known,  however,  this  is  not  necessary  if  the  effective 


2-23 


channel  order  L  is  known,  which  is  also  a  known  constant  in  the  above  subspace 
channel  identification  derivation  [4:697].  Accurate  channel  length  estimation  is  also 
critical  to  extracting  the  equalizer  tap  coefficients  from  the  channel  estimate,  with  a 
poor  estimate  leading  to  an  unstable  solution.  Channel  order  estimation  is  preferred 
over  noise  power  estimation  since  numerical  methods  exist  for  estimating  the  effective 
channel  order  without  knowledge  of  the  noise  variance  [8].  Finally,  like  CMA,  the 
subspace  algorithm  relies  on  the  linear  independence  of  the  virtual  channels  [20]. 

2-4-6  Spectrum  Fitting  Algorithms.  Spectrum  fitting  channel  identification 
algorithms  based  on  second  order  statistics  are  a  frequency  domain  parallel  to  the 
time  domain  block  subspace  algorithms.  As  an  SOS-based  algorithm,  spectrum  fitting 
has  the  potential  to  accurately  identify  the  channel  filtering  coefficients  in  relatively 
few  symbol  periods  as  compared  to  HOS-based  algorithms.  Like  the  block  subspace 
method,  spectrum  fitting  performs  matrix  operations  on  a  collected  set  of  fractionally 
sampled  symbols,  making  it  periodically  adaptive.  The  channel  properties  that  ensure 
a  unique  identification  solution  are  the  same  as  those  for  the  block  subspace  algorithm. 
The  frequency  domain  approach  provides  an  alternative  formulation  of  the  conditions 
under  which  SOS-based  channel  identification  succeeds  [3,19].  Of  note,  Ding  and  Li 
point  out  that  this  method  is  is  not  affected  by  the  kurtosis  of  the  input  signal,  unlike 
most  CMA  algorithms  which  break  down  for  near-Gaussian  sources  [3]. 

The  general  strategy  of  spectrum  fitting  algorithms  relies  on  an  estimate  of  the 
power  spectral  density  of  the  received  signal.  In  the  multichannel  model,  the  data 
sequence  in  each  branch  (or  antenna)  is  stationary,  so  autocorrelation  can  be  estimated 
via  time  averaging.  Taking  the  Fourier  transform  of  each  antenna  produces  O  power 
spectral  density  functions,  referred  to  as  “cyclic  spectra”  [19].  The  coefficients  of  the 
cyclic  spectra  can  be  used  to  construct  a  block  convolution  matrix  similar  in  spirit  to  G 
in  the  subspace  decomposition  algorithm.  This  leads  to  an  error  minimization  solution 
involving  a  singular  value  decomposition  of  the  block  convolution  matrix.  Because 
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spectrum  fitting  is  a  frequency  domain  analog  of  block  subspace  decomposition,  it  is 
not  investigated  further  in  this  thesis. 

2-4-7  Cumulants  Method.  Shalvi  and  Weinstein  have  proposed  a  direct  blind 
equalization  method  based  on  matching  selected  second  and  fourth  order  statistics  of 
the  channel  input  and  equalization  filter  output.  The  result  is  an  adaptive  gradi¬ 
ent  descent  algorithm  similar  in  spirit  to  CMA,  with  the  advantage  that  source  signal 
distributions  with  kurtosis  exceeding  the  kurtosis  of  the  Gaussian  distribution  (super- 
Gaussian)  are  allowable.  Specifically,  the  Shalvi- Weinstein  method  assumes  the  source 
distribution  consists  of  zero-mean,  identically  distributed  independent  symbols  and 
that  the  first  through  fourth  moments  of  the  source  distribution  exist.  Additionally, 
the  excess  kurtosis  of  the  source  distribution  cannot  equal  zero,  so  Gaussian  distri¬ 
butions  are  not  equalized  by  this  method.  The  channel  constraints  for  the  Shalvi- 
Weinstein  algorithm  are  the  same  as  those  for  CMA  [14]. 

The  mathematical  foundation  developed  in  [14]  upon  which  this  algorithm  is 
based  is  the  following  theorem: 

If  £[|i(n)|2]  =E[\x(n)\2}  , 

then  \K  (x(n))\  <  \K  (x(n))\  ,  (2.40) 

where 

K(x{n))  =  E[\x(n)\*]  -  E2  [\x(n)\2]  -\E[x2(n)]\2  .  (2.41) 

K  (x(n))  is  kurtosis  as  defined  in  [14].  The  kurtosis  magnitude  of  the  input  and 
output  symbols  in  (2.40)  are  equal  when  the  combined  channel/filter  response  is  a 
discrete  unit  impulse.  In  other  words,  when  the  output  sequence  of  the  equalizing 
filter  is  identical  to  the  input  sequence  with  some  delay,  then  the  kurtosis  of  the  two 
sequences  is  equal  (assuming  the  sequences  have  equal  power)  [14], 

This  leads  to  the  adaptive  equalization  strategy  of  selecting  filter  coefficients 
to  maximize  the  kurtosis  of  the  output  sequence  under  the  equal  power  constraint 
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E  [|a?(rz) | 2]  =  E  [|x(n)|2] .  The  constrained  maximization  problem  leads  to  a  gradient 
ascent  formulation.  This  gradient  ascent  form,  in  turn,  can  be  transformed  into  a 
symbol-by-symbol  iterative  approximation: 

w(n  +  1)  =  w (n)  +  /i  ■  [sgn  K(x(n))\  \x(n)\~  x(n)y*(n)  .  (2.42) 

where  y *(n)  is  the  received  symbol  stream  passed  through  a  whitening  filter.  For 
source  signals  with  sub-Gaussian  kurtosis,  the  tap  update  algorithm  proposed  in  [14] 
simplifies  to  gradient  descent  using  the  CMA  error  term  (2.29).  Because  the  DPSK 
constellations  used  in  this  thesis  are  sub-Gaussian,  the  cumulants  method  is  not  stud¬ 
ied  as  a  potential  non-cooperative  receiver. 
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III.  LPI  Signal  Design 


he  goal  of  this  thesis  is  to  investigate  a  communications  signal  structure  that 


_L  resists  equalization  by  blind  methods,  leveraging  non-cooperative  receivers’  re¬ 
liance  on  blind  equalization  to  create  an  LPI  signal.  In  designing  such  a  scheme,  the 
impact  to  the  cooperative  communications  link  must  also  be  considered,  both  in  terms 
of  transmission  accuracy  and  system  complexity.  This  chapter  proposes  a  strategy  of 
introducing  a  time  dither  into  the  symbol  synchronization  of  the  transmitter.  This 
dither  creates  dynamic  ISI  at  the  receiver  input  which  hinders  blind  equalization.  A 
cooperative  receiver,  however,  can  successfully  deconvolve  the  symbol  sequence  if  the 
timing  dither  of  the  symbol  transmissions  is  known.  The  signaling  scheme  is  developed 
in  detail  in  Section  I5A1  Section  [5121  presents  cooperative  receiver  designs  for  this  com¬ 
munications  signal  and  Section  15131  discusses  the  behavior  of  notional  non-cooperative 
receivers  with  blind  algorithms  with  respect  to  the  symbol  dither  signal. 

3.1  Signal  Design 

All  equalization  techniques  discussed  in  Chapter  II  rely  on  a  sufficiently  slow 
time  variation  in  the  effective  channel  impulse  response.  This  thesis  investigates  the 
feasibility  of  achieving  LPI  characteristics  in  a  communications  signal  by  injecting 
time-varying  ISI  into  the  waveform  via  dither  in  the  symbol  transmission  synchro¬ 
nization.  To  prevent  signal  interception,  the  resulting  ISI  must  be  severe  enough 
to  prevent  meaningful  demodulation  and  vary  in  time  in  a  manner  that  prohibits 
effective  blind  equalization. 

More  explicit  discussion  of  synchronization  dither  requires  additional  notation. 
In  most  traditional  transmitters,  communications  symbols  are  transmitted  at  a  con¬ 
stant  rate  with  the  leading  edge  of  each  symbol  exactly  Tsym  seconds  apart.  To  create 
ISI  in  the  LPI  waveform,  let  the  leading  edge  of  each  symbol  be  delayed  up  to  an  entire 
additional  symbol  period.  In  other  words,  the  LPI  waveform  is  composed  of  a  series 
of  time  blocks  of  duration  Tsym ,  each  block  containing  the  leading  edge  of  exactly 
one  communications  symbol  somewhere  within  the  block.  This  dithering  method  pre- 
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x[n) 


x^ct\t) 


Figure  3.1:  A  dither  block  in  the  transmitter  creates  ISI. 


serves  symbol  order  and  fixes  the  average  symbol  transmission  rate  to  a  constant  rate 
of  1  /Tsym,  simplifying  the  cooperative  receiver  demodulation  process. 

The  delay  of  the  nth  symbol  relative  to  the  leading  edge  of  the  nth  symbol  period 
is  d(n).  This  delay  is  applied  after  the  DAC  in  the  single  channel  communications 
model  in  Figure  12711  as  shown  in  Figure  13.11  This  delay  implementation  of  dither 
requires  a  new  expression  for  x^ct\t)\ 

TV— 1 

X^\t)  =  ^2  x(n)S  (t  -  nTsym  -  d(n))  ,  (3.1) 

71=0 

where  d{n )  is  a  sample  of  random  variable  D  with  a  probability  density  function 
p2\t)  and  ct  denotes  a  continuous  time  waveform.  This  thesis  considers  uniform 
density  functions,  although  more  exotic  density  functions  are  possible,  so  long  as  the 
function  is  zero  for  t  <  0  and  t  >  Tsym.  The  uniform  dither  distribution  is  completely 
defined  by  the  dither  range,  which  is  0  <  d(n)  <  Tmax.  The  probability  density 
function  can  be  written  in  terms  of  Tmax  according  to 


p'S'W 


0<t<  Tm„ 

0,  t  <  0  or  t  >  Tmax  . 


(3.2) 


Clearly,  for  p^\t)  to  satisfy  the  restrictions  above,  the  dither  range  Tmax  is  restricted 
to  [0, Tsym\. 

In  practice,  the  sequence  of  symbol  delays  is  a  pseudorandom  sequence  with 
terms  selected  from  a  discrete  set  of  possible  delays.  In  this  discrete  case,  the  symbol 
period  is  subdivided  into  P  evenly-spaced  possible  dither  values.  The  probability 
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density  function  for  a  discrete  uniform  random  variable  D  is  given  by 


(3.3) 


where  pmax  must  be  from  the  set  of  integers  ranging  over  [0,  P ).  The  parameter  pmax 
in  the  discrete  distribution  is  related  to  the  dither  range  Tmax  according  to 


Pmax 


PT 

l  ±  ry 


Tc 


sym 


(3.4) 


where  [xj  is  the  largest  integer  not  exceeding  x. 

As  a  pseudorandom  sequence,  the  delay  code  d{n)  will  have  some  finite  pe¬ 
riod  Nd.  In  most  cases,  this  thesis  considers  long  delay  codes  with  periods  Nd  larger 
than  the  number  of  transmitted  symbols.  The  cooperative  receiver  described  in  Sec¬ 
tion  13.2.11  however,  requires  short  period  relative  to  the  number  of  symbols  available 
for  training  the  equalization  hlter(s). 

The  severity  of  ISI  in  this  waveform  is  a  function  of  the  pulse  shape  and  the 
symbol  delay  probability  density  function.  Conceivably,  symbol-to-symbol  variation 
in  the  pulse  shape  of  the  communications  waveform  could  be  used  as  an  alternative 
or  in  addition  to  symbol  synchronization  dither  to  inject  ISI.  Pulse  shape  variation, 
however,  introduces  additional  complexity  on  matched  filtering  in  the  cooperative  re¬ 
ceiver.  For  this  reason,  the  LPI  benefits  due  to  intentional  ISI  are  characterized  using 
a  fixed  pulse  shape  and  varying  symbol  dither  according  a  pseudorandom  sequence. 


3.2  Cooperative  Receiver  Strategies 

This  thesis  develops  cooperative  receiver  designs  by  applying  signal  processing 
techniques  on  the  sampled  matched  filter  output.  To  effectively  demodulate  a  pseudo¬ 
random  dither  waveform  (PDW),  a  cooperative  receiver  uses  some  knowledge  of  the 
delay  sequence  used  to  generate  the  waveform  in  the  transmitter.  Two  methods  are 
suggested  here.  In  the  first,  the  periodic  nature  of  the  symbol  delay  sequence  is  uti- 
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Figure  3.2:  Generic  structure  for  the  cooperative  and  non-cooperative  receivers 

described  in  this  chapter. 


lized.  In  the  second  cooperative  receiver  design,  full  knowledge  of  the  delay  sequence 
is  assumed  to  construct  an  equalization  filter  according  to  the  MSE  algorithm  (in  a 
line-of-sight  environment). 

Both  the  cooperative  and  non-cooperative  receivers  discussed  in  the  remainder 
of  this  chapter  use  the  structure  shown  in  Figure  3.21  The  signal  at  the  receiver 
input  r^ct\t)  is  passed  through  a  matched  filter  based  on  the  communication  signal 
pulse  shape.  For  non-cooperative  receivers,  it  is  assumed  that  the  pulse  shape  can 
be  estimated  using  the  power  spectral  density  of  the  received  waveform.  The  rate 
at  which  the  matched  filter  output  is  sampled  and  the  Post-Processing  block  are 
specific  to  the  different  receiver  structures.  The  resulting  estimate  of  the  source 
symbol  sequence  x  (n)  is  demodulated  to  produce  the  output  bit  stream. 


3.2.1  Bank  of  Equalizers.  This  equalization  structure  takes  advantage  of  the 
periodicity  of  the  symbol  delays  in  the  delay  code  sequence.  Assuming  the  physical 
channel  is  static,  all  samples  at  the  matched  filter  ouptut  spaced  ( NdTsym )  seconds 
apart  are  subject  to  the  same  effective  channel.  Thus  each  collection  of  these  samples 
is  WSS.  A  cooperative  receiver  with  knowledge  of  the  delay  sequence  period  can  cycle 
through  a  bank  Nd  equalization  filters,  one  for  each  effective  channel.  Thus,  this 
multiple  filter  approach  is  only  feasible  for  dither  codes  with  small  periods  due  to 
convergence  and  coefficient  storage  considerations.  The  impact  of  short  period  codes 
is  discussed  in  greater  detail  below.  The  receiver  structure  for  a  fractionally  spaced 
equalizer  with  oversampling  rate  2  and  delay  code  period  Nd  is  depicted  in  Figure  13131 
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Figure  3.3:  By  applying  a  different  equalizer  for  each  pseudorandom  dither 

entry  in  the  dither  code,  the  effective  channel  at  each  equalizer  input  is  static,  allowing 
the  filters  to  converge. 


The  remaining  task  is  to  find  the  equalizer  filter  taps  such  that  the  output  of  each 
filter  is  an  estimate  of  a  unique  decimation  of  the  original  symbol  sequence  x(n  —  5).  In 
other  words,  to  produce  a  complete  estimate  x(n),  each  filter  must  estimate  x(n  —  5) 
decimated  by  Nj,  but  with  different  phases.  In  addition,  the  phase  of  each  must  be 
known  so  that  filter  outputs  can  be  recombined  in  the  correct  order.  Because  the 
number  of  (fractionally  spaced)  equalizer  filters  matches  the  dither  sequence  period, 
the  sample  statistics  at  each  filter  input  are  stationary  and  many  classical  equalization 
algorithms  are  applicable. 

Of  the  equalization  algorithms  described  in  Chapter  HH  the  trained  approach  of 
adaptive  LMS  is  the  best  candidate.  The  adaptive  LMS  algorithm  can  be  directly 
applied  in  each  filter  branch  to  train  each  equalizer  without  introducing  any  increase 
in  algorithm  complexity.  The  simplest  method  to  ensure  that  each  estimated  symbol 
sequence  at  the  equalizer  outputs  represent  unique  subsets  of  the  transmitted  symbol 
sequence  x(n)  is  to  use  the  same  5  parameter  for  each  filter.  With  this  approach,  the 
demultiplexer  can  recombine  the  filter  outputs  using  the  same  order  of  parsing  used 
in  the  multiplexer. 
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The  CMA  method  does  not  extrapolate  to  a  multiplexed  approach  quite  as  eas¬ 
ily.  The  difficulty  arises  from  the  indeterminate  nature  of  the  delay  6  between  the 
source  sequence  and  the  estimated  sequence  at  the  equalizer  output  in  the  classical 
CMA  implementation.  This  algorithm  approximates  the  sequence  x{n  —  5),  but  5  is 
not  specified  explicitly.  In  the  equalizer  bank  approach  proposed  here,  each  equal¬ 
izer  is  computed  independently  and  may  result  in  different  integer  values  of  8.  If 
this  happens,  the  parallel  equalizers  may  produce  estimates  of  the  same  decimation 
of  the  original  symbol  sequence,  while  other  subsequences  are  not  estimated  at  all. 
Notionally,  it  should  be  possible  to  force  each  equalizer  to  produce  a  symbol  sequence 
uncorrelated  with  the  output  of  each  other  equalizer,  but  this  requires  increased  com¬ 
plexity.  Furthermore,  even  if  the  symbol  subsequence  in  each  branch  is  independent, 
the  proper  order  for  recombining  the  equalizer  outputs  into  a  single  sequence  is  not 
known. 

The  subspace  algorithm  first  computes  a  channel  estimate  which  is  used  to  find 
an  equalization  filter,  so  the  relative  delay  8  between  the  input  and  output  symbol 
sequences  can  be  controlled  explicitly.  The  robustness  of  the  equalizer  with  respect  to 
noise,  however,  is  highly  dependent  on  the  value  of  8,  so  only  a  subset  of  possible  delays 
are  reasonable  candidates  for  a  particular  effective  channel  impulse  response  [11].  It 
is  conceivable  that  the  effective  channels  seen  at  the  equalizer  inputs  will  be  similar 
enough  to  have  a  common  6  value  that  is  near-optimal  for  each  filter  branch,  although 
not  guaranteed. 

In  addition  to  only  minor  additional  complexity  in  the  receiver  architecture,  ad¬ 
vantages  of  the  LMS  bank  of  equalizers  scheme  include  intrinsic  multipath  mitigation 
and  modest  a-priori  knowledge  of  the  delay  sequence  in  the  transmitter.  The  only 
knowledge  of  the  pseudorandom  delay  sequence  used  to  create  the  waveform  in  the 
transmitter  is  the  period  Nj.  Any  change  in  the  delay  code  or  the  physical  channel 
merely  requires  the  equalization  filters  to  be  retrained.  Typically,  a  communications 
system  using  LMS  equalization  will  retransmit  the  training  sequence  periodically.  If 
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the  sequence  is  sufficiently  long  for  the  filters  in  the  receiver  to  converge,  the  dither 
code  could  be  changed  in  the  transmitter  every  time  the  training  sequence  is  resent. 

A  major  limitation  of  this  method  is  slow  filter  convergence  and  increased  mem¬ 
ory  requirements  to  store  coefficients  for  Nd  filters.  The  number  of  training  symbols 
for  convergence  is  increased  by  a  factor  of  Nd  over  a  standard  transmitter /receiver 
configuration  because  each  equalizer  must  be  trained  independently  and  only  sees 
one  symbol  for  every  Nd  symbols  transmitted.  For  long  delay  codes  (with  large  pe¬ 
riods),  equalization  filter  convergence  time  may  be  prohibitively  large.  Similarly,  the 
coefficients  for  each  of  the  Nd  filters  must  be  stored  within  the  receiver,  limiting  the 
code  period.  Certainly,  any  decrease  in  the  dither  code  length  may  facilitate  code 
estimation  by  a  non-cooperative  receiver.  In  addition,  a  short  code  uses  few  samples 
of  the  random  variable  D,  and  may  not  be  a  good  representation  of  the  governing  dis¬ 
tribution.  As  a  result,  a  particular  (short  period)  code  may  have  unique  ISI  severity 
corresponding  to  widely  disparate  degrees  of  LPI  protection. 

3.2.2  MSE  Minimization  Using  Known  Delay  Sequence.  If  the  pulse  shape 
and  the  delay  code  sequence  of  the  PDW  are  known  by  the  receiver,  the  ISI  at  the 
matched  filter  output  due  to  overlapping  symbol  shapes  can  be  expressed  as  a  known 
time- varying  impulse  response.  In  a  line-of-sight,  lossless  environment,  each  sample  of 
the  matched  filter  output  can  be  written  in  terms  of  a  discrete  impulse  response  that 
describes  the  contribution  of  input  symbols  to  the  output  sample  in  question.  For  a 
completely  known  channel,  MSE  minimization  is  a  common  approach  for  computing  a 
discrete  FIR  filter  to  deconvolve  the  channel  effects.  Typically,  the  MSE  minimization 
algorithm  assumes  a  constant  channel  impulse  response  and  produces  a  single  filter. 
In  the  PDW  equalization  problem,  a  unique  MSE  minimization  solution  must  be 
found  for  each  of  the  Nd  elements  in  the  delay  code  sequence. 

The  theoretical  foundation  of  the  MSE  algorithm  presented  in  [18:287]  is  sim¬ 
ilar  to  the  adaptive  LMS  development  discussed  in  Section  12.4.21  Instead  of  using 
stochastic  methods  to  estimate  the  expected  value  of  the  correlation  matrices  in  the 
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cost  function  (]3.5lh  however,  the  matrix  coefficients  of  Ry  and  pxy  are  known. 


J(n)  =  E  [x*[n  —  S)x(n  —  5)]  —  2Re  {e  \x(n  —  S)yH(n )]  w*} 

+  wrE  [y(n)yH(n)j  w* 

=  E  [x*(n  —  8)x(n  —  5)]  —  2Re  {p^w*}  +  wT Ryw*  (3.5) 

By  setting  the  cost  function  gradient  (3.6)  to  zero,  the  optimal  equalization  filter  in 
the  mean-square  error  sense  can  be  found  explicitly  (377) . 


VwJ(u)  =  -2p^y  +  2wrRy 


(3.6) 


W  opt  = 


(3.7) 


In  adapting  the  MSE  algorithm  for  the  pseudorandom  delay  waveform,  a  frac¬ 
tionally  spaced  implementation  is  considered.  The  vectorized  matched  filter  output 
y(n)  is  a  columnwise  stack  of  the  sampled  subchannels  at  the  matched  filter  output 
(see  (2.17)  and  (2.22)  in  Section  [2.3.21).  Similarly,  the  subchannel  equalization  filters 
are  stacked  in  a  column: 

w  (°)(n) 


w  (n)  = 


(3.8) 


In  traditional  MSE,  the  effective  channel  impulse  response  is  static  and  the 
source  sequence  x(n)  is  uncorrelated,  so  the  received  sequence  within  each  subchannel 
y^\n)  is  stationary  and  the  autocorrelation  E  [y^ (ni)y^* (n2)]  is  dependent  only  on 
the  difference  5  =  ri2  —  n\  and  reduces  to  E  [y^\n)y^* (n  —  5)] .  In  the  pseudorandom 
delay  waveform,  ISI  due  to  pulse  shape  dither  creates  different  effective  channels. 
To  signify  this,  let  H^(n)  denote  the  convolution  matrix  that  maps  the  transmitted 
symbols  in  vector  x(n)  to  the  output  sample  vector  y(n).  Note  that  the  convolution 
matrix  must  now  be  indexed  by  n  to  account  for  the  potentially  unique  delay  applied 
to  each  symbol.  The  autocorrelation  matrix,  then,  is  also  dependent  on  n  according 
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to 


RyW  =  E  [y(«)y H{n)]  , 


=  E 


(Hjv(n)x(n)  +  b(n))  (Hjv(n)x(n)  +  b(n))H 
=  E  [HAr(n)x(n)x^(77,)H^(n)]  +  E  [b(n)b^(n)]  , 
=  Hyv(n)R;(.H^(n)  +  Rb  . 


(3.9) 


In  a  similar  manner,  the  cross-correlation  vector  p xy(n)  is  dependent  on  the  symbol 
index  n: 


pxy(n)  =  E[y*(n)x(n- 5))  , 

=  E  [(Hjv(n)x(n)  +  b(n))*  x[n  -  5)]  , 

=  H *N(n)E  [x*(n)x(n  —  5)}  .  (3.10) 

Finally,  the  MSE  optimal  filter  is  now  indexed  by  symbol  as  well: 

wopt{n)  =  (Ry(n))  1pxy(n).  (3.11) 

The  remaining  task  is  to  find  the  parameters  that  define  Ry(n)  and  p XJ/(n). 
Specifically,  RN(n),  R^,,  R;,,  and  E  [x.*(n)x(n  —  5)]  must  all  be  assigned  numeric 
values. 

If  the  source  sequence  x(ri)  is  zero- mean  iid,  then  Ra  is  simply  the  identity 
matrix  scaled  by  E  [x*(n)x(n)]  and  E  [x*(n)x(n  —  5)]  has  only  one  non- zero  term 
with  value  E  [x*(n)x(n)]  at  the  (6  +  l)th  element. 

The  elements  within  noise  autocorrelation  matrix  Rb  can  be  expressed  analyt¬ 
ically  in  terms  of  the  white  Gaussian  noise  power  spectral  density  in  the  physical 
channel  N0  and  the  matched  filter  fit).  Each  element  of 

R b  =  E  [h{n)bH(n)]  (3.12) 
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can  be  expressed  as  E  [b'(ki)b'*(k2)]-  Substituting  for  the  definition  of  b'(k)  from  (2.7) 
in  Section  12.3.11  gives 


E[b'(k1)b'*(k2)]  = 

JZ,  f{a)  (n)  v{a)  (T,  -  n)  dn  JZ  f(cth  (75)  (r2  -  T2)  <iT2 


(3.13) 


where 


T\  = 

r2  = 


fciT. 


sym 


o 

k-2  Tsym 


o 


The  whiteness  property  of  the  AWGN  signal  u^ct\t)  implies  that  the  autocorre¬ 
lation  is 

E  =  {  N°'  ^  ~  h  (3.14) 

[0,  ti  =£  t2  ■ 

Applying  this  to  (3. 13]).  the  two-dimensional  integration  over  area  dTidr2  is  only  non¬ 
zero  when  Ti  —  Ti  =  T2  —  r2.  As  a  result,  (13.13)  reduces  to  a  line-integral  of  the 
form 


/OO 

/(^(rO/^^-Ti  +  rOdn  .  (3.15) 

-OO 

The  integrand  of  (13.15)  is  simply  the  convolution  of  pulse  shape  f^ct\t)  and  the 
matched  filter  /(rf)*(_ t).  Thus,  the  elements  of  Rj,  are  samples  of  the  convolu¬ 
tion  (13.16).  scaled  by  the  noise  power  spectral  density  N0  (13.17). 


a{a\t)  =  *  /<rfl'(-«) 


f(Ct>  (T)fW(t  +  T)dT 


p(ci)* 


(3.16) 
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E  [b\k1)b,’(k2)]  =  Nja)  [(n2  -  m)Tnm  +  (i,  -  ii)T,ym/0] 


(3.17) 


As  an  alternative  to  analytically  solving  for  the  elements  of  R&  as  shown  above, 
the  filtered  noise  correlations  can  be  estimated  by  time-averaging  the  samples  of 
the  matched  filter  output  when  the  communications  waveform  is  not  present,  taking 
advantage  of  the  WSS  nature  of  the  filtered  noise. 

Finally,  the  convolution  matrix  Hat(t7,)  must  be  found.  As  mentioned  above, 
H/v(n)  is  indexed  by  symbol  n  because  of  the  unique  timing  delay  of  each  symbol.  In 
a  convolution  matrix,  each  column  represents  the  sampled  lattice  of  the  continuous 
time  effective  channel  impulse  response  that  models  the  amplitude  scale  factor  of  a 
particular  symbol  at  future  points  in  time.  In  the  PDW  model  described  in  this  thesis, 
the  continuous  time  effective  impulse  response  (t)  is  the  same  for  each  transmit¬ 
ted  symbol.  The  pseudorandom  timing  delay  of  each  symbol,  however,  impacts  the 
synchronization  of  the  sampling  lattice  relative  to  h^jj  (t)  that  models  the  influence 
of  a  symbol  at  each  sample  of  the  matched  filter  output.  In  other  words,  each  column 
of  H ;v(n)  contains  the  a  sampled  lattice  of  the  effective  channel  impulse  response, 
but  with  potentially  differing  sampling  phases  as  determined  by  the  delay  of  each 
particular  symbol.  Taking  the  sampling  lattice  phase  offset  due  to  a  symbol’s  delay 
d(n),  the  Ah  subchannel  impulse  response  defined  in  (}2.14j)  is  now  written  as 

hw(n)  =  [h{eff  (( i  -  1  )Tsym/0  -  t0  -  d(n )) ,  h(eff  ((O  +  i  -  l)Tsym/0  - 10  -  d(n )) , 
h^f  ((20  +  *  -  1  )Tsym/0  -  t0  -  d(n)) ,  ..., 
hg}  (((L  -l)0  +  i-  1  )Tsym/0  -  t0  -  d(n))]T  , 

=  [h®( 0),  h®(  1),  2),  ...  h®{L  -  1)]T  .  (3.18) 
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Consequently,  the  subchannel  convolution  matrix  for  the  pseudorandom  delay 
waveform  is 

H$(n)  = 

N+L 

/ - ^ - \ 

"^(0)  •••  h®L+1(L-  1)  0  .  0 

0  0)  •••  h®L(L- 1)  0  0 

N 

0  ■"  0  ^n-N+ 1(0)  ■"  ~  1) 

(3.19) 

and  the  composite  convolution  matrix  Hjv(n)  is  constructed  by  stacking  the  subchan¬ 
nel  convolution  matrices  column-wise  as  in  (2.22).  With  these  definitions,  w opt(n)  can 
now  be  calculated  from  (3.11). 

Because  the  equalization  filters  are  deterministic,  this  modified  MSE  cooperative 
receiver  design  has  no  need  to  estimate  signal  statistics  like  the  algorithms  presented 
in  Chapter  [III  Instead,  the  equalization  filter  taps  are  computed  analytically.  For 
a  static  physical  channel  and  a  particular  delay  code  sequence,  Nd  filters  must  be 
calculated.  Thus,  the  initial  computational  complexity  for  a  new  effective  channel 
involves  computing  Nd  matrix  inverses,  which  only  need  to  be  computed  once. 

The  symbol  estimation  of  the  modified  MSE  receiver  is  optimal  in  the  mean- 
squared  error  sense  with  respect  to  the  available  samples.  Because  the  matched 
filter  output  is  sampled  at  a  constant  rate  and  samples  are  not  synchronized  to  the 
symbol  dithers,  some  symbol  energy  loss  is  expected  due  to  non-optimal  sampling. 
Oversampling  helps  compensate  for  this  suboptimal  sampling  strategy. 

Computation  of  the  effective  channel  impulse  response  is  not  directly  addressed 
in  the  receiver  design.  Under  line-of-sight  conditions,  knowledge  of  the  pulse  shape 
(which  is  assumed  to  allow  the  receiver  to  use  a  matched  hlter)  and  the  transmission 
delay  is  sufficient.  In  the  case  of  a  multipath  physical  channel,  the  channel  model 
must  be  estimated,  possibly  via  transmission  of  a  known  training  waveform. 
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{  bout } 


Figure  3.4:  Generic  structure  for  the  cooperative  and  non-cooperative  receivers 

described  in  this  chapter. 


The  calculation  of  each  equalization  filter  w opf(n)  requires  that  the  cooperative 
receiver  have  full  knowledge  of  the  pseudorandom  symbol  delay  sequence  used  in  the 
transmitter.  Also,  the  receiver  must  be  fully  synchronous  in  the  sense  that  the  symbol 
estimate  x(n  —  5)  relics  on  the  application  of  the  corresponding  filter  w opt(n). 


3.3  Non-cooperative  Receivers 

Several  non-cooperative  receiver  structures  are  considered  here  with  respect  to 
their  ability  to  exploit  an  intercepted  PDW.  Each  notional  non-cooperative  receiver 
is  assumed  to  have  knowledge  of  the  symbol  period,  pulse  shape,  and  symbol  con¬ 
stellation,  but  no  knowledge  of  the  symbol  dither  constants  d(n).  A  basic  receiver 
that  makes  no  attempt  to  account  for  the  symbol  timing  dither  is  be  considered  first. 
Next,  a  receiver  with  blind  multipath  equalization  processing  is  considered.  A  third 
possible  non-cooperative  receiver  strategy  involves  estimating  the  symbol  delays  to 
sample  the  matched  filter  output  asynchronously. 

3.3.1  Basic  Receiver.  This  section  considers  the  receiver  shown  in  Fig¬ 
ure  13.4"  Note  that  the  matched  filter  output  is  sampled  at  a  rate  of  1/Tsym,  with  con¬ 
stant  time  Tsym  between  successive  samples.  No  post-processing  is  performed.  The 
sampling  phase  is  synchronized  to  the  optimum  time  for  a  symbol  with  delay  d(n) 
equal  to  the  mean  of  the  probability  density  function  associated  with  D. 

Assuming  a  line-of-sight  physical  channel  and  a  traditionally  modulated  wave¬ 
form  (with  no  symbol  timing  dither),  if  the  matched  filter  output  samples  are  syn¬ 
chronized  to  the  received  waveform,  intersymbol  interference  is  not  present.  For  the 
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PDW  in  a  line-of-sight  channel,  the  ISI  is  now  solely  due  to  pulse  shape  overlap  caused 
by  symbol  dither.  Thus,  any  statistically  significant  deviation  from  the  bit  error  rate 
without  symbol  timing  dither  can  be  attributed  to  the  ISI  induced  by  the  PDW. 

3.3.2  Blind  Equalizer  Receiver.  Because  the  structure  of  multipath  inter¬ 
ference  is  similar  to  the  ISI  caused  by  symbol  timing  dither  in  the  PDW,  a  non- 
cooperative  receiver  using  blind  equalization  techniques  is  considered  here.  Equal¬ 
ization  to  mitigate  the  ISI  within  the  PDW  is  attempted  using  equalizer  filters  de¬ 
termined  by  the  Constant  Modulus  Algorithm  and  the  subspace  decomposition  algo¬ 
rithm.  Fractionally-spaced  versions  of  both  equalization  algorithms  are  implemented, 
with  the  matched  filter  output  in  Figure  fOI  sampled  at  a  rate  of  2/Tsym.  As  in  the 
basic  receiver  architecture  described  in  the  previous  section,  the  sampling  is  synchro¬ 
nized  to  the  mean  symbol  dither. 

Equalization  using  CMA  and  subspace  decomposition  is  described  in  Sectionl2.4.3l 
and  Section  I2.4.51  respectively.  Both  of  these  methods  rely  on  the  cyclostationarity 
of  the  oversampled  matched  filter  output,  which  is  clearly  violated  by  the  dynamic 
ISI  within  the  PDW.  However,  it  is  conceivable  that  the  average  interference  between 
symbols  based  on  their  relative  indexes  could  lead  to  some  benefit  from  an  equaliza¬ 
tion  filter.  Any  potential  BER  improvement  from  such  a  strategy  would  depend  on 
the  particular  pulse  shape  and  the  distribution  of  the  delay  terms  d{n). 

3. 3. 3  Symbol  Dither  Estimation  Receiver.  In  the  previous  two  non-cooperative 
receivers,  standard  techniques  are  considered  with  respect  to  the  PDW.  The  archi¬ 
tecture  considered  here  specifically  targets  interception  of  a  PDW  (in  a  line-of-sight 
environment  or  known  physical  channel  impulse  response).  Assuming  knowledge  of 
the  pulse  shape,  signal  power,  and  the  symbol  constellation,  a  symbol  can  be  assigned 
to  each  symbol  period  that  best  fits  the  received  waveform  in  a  sum  (or  mean)  squared 
error  sense.  Within  each  symbol  period,  the  receiver  searches  for  the  symbol  and  de¬ 
lay  pair  that  minimizes  the  sum  squared  difference  between  the  received  waveform 
and  the  estimated  symbols.  This  is  accomplished  by  creating  an  estimated  waveform 
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in  the  receiver  via  modulation  of  the  estimated  symbol  sequence  and  application  of 
the  estimated  time  dither  to  each  symbol. 

Symbol  and  delay  estimation  proceeds  as  follows.  Using  the  known  effective 
channel  impulse  response  and  signal  power  at  the  receiver  input,  the  continuous  time 
signal  representing  each  symbol  in  the  constellation  can  be  created  within  the  receiver. 
Using  complex  baseband  symbols,  let 

■f’M  =  (3.20) 


be  the  continuous  time  shape  for  complex  symbol  s*  at  the  receiver  (see  (14.2|)).  The 
receiver  must  assume  some  finite  set  of  Q  possible  symbol  dither  values  for  each  of 
N  symbols.  The  initial  estimated  waveform  is  y^l\t )  =  0.  Symbols  are  estimated  by 
taking  N  iterations  of  the  three  steps  below. 


(1)  :  z)  (t,  i,  q,  n )  =  y)_[  ( t )  +  s\  t-  qTsym/Q  -  nTt 


\w i  "i  'll  —  '  “i  \v  syra /  ^  •v-i-sym 

(2)  :  ij,  qjUj  =  arg min  <j  j^+NT^  yW  (t)  -  zft]  (t,  i,  q,  n) 

i,q,n  f  L 

(3)  :  yft]  (t)  =  yW  (t)  +  s[f  (t  -  qfTsym/Q  -  nfTsym 


dt 


(3.21) 


In  step  (1),  all  valid  combinations  of  symbol  and  symbol  are  used  to  create  temporary 
signal  estimates  by  adding  the  potential  symbols  to  the  previous  iteration’s  estimate. 
Step  (2)  selects  the  symbol  i,  symbol  index  n,  and  symbol  dither  q  set  that  results  in 
the  smallest  sum  error  term.  Finally,  the  new  signal  estimate  is  created  in  step  (3). 


The  ranges  of  i  and  q  are  the  same  for  each  iteration.  All  symbols  in  the  M— ary 
constellation  are  considered  for  i  and  the  dither  estimate  ranges  over  the  set  of  all 


possible  integer  steps:  q  G 


0...Q-1 


The  search  range  of  n  begins  as  the  integers 
from  0  to  N  —  1,  inclusive,  but  shrinks  with  each  iteration.  After  each  iteration, 
one  symbol  is  estimated  and  no  future  iterations  can  assign  a  symbol  to  that  symbol 
period,  so  rij  is  removed  from  the  search  range.  After  N  iterations,  one  symbol 
has  been  estimated  per  each  symbol  period.  Note  that  this  process  assumes  that 
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the  receiver  is  able  to  synchronize  to  the  symbol  period  boundaries  of  the  received 
waveform. 


The  computational  complexity  of  this  approach  is  not  trivial.  The  search  space 
of  the  first  iteration  alone  requires  computing  ^ MQN^j  error  terms,  where  N  will  often 
be  quite  large.  To  be  feasible,  only  a  fraction  of  the  received  symbols  is  considered  at 
a  time. 


Optimization  of  this  symbol  estimation  technique  is  beyond  the  scope  of  this 
thesis.  The  iterative  method  described  here  could  be  replaced  by  a  Viterbi-like  tech¬ 
nique  the  considers  (^MQ'j  branches  at  each  node.  Sparse  channel  identification 
techniques  may  also  apply  to  this  problem  by  considering  the  series  of  symbols  (mod¬ 
eled  as  complex  impulses  in  the  PDW)  as  the  impulse  response  of  a  sparse  multipath 
channel.  Symbol  estimation  is  then  the  equivalent  to  the  sparse  channel  identifica¬ 
tion  problem,  with  the  additional  restrictions  that  each  channel  ray  must  be  one  of 
the  complex  symbols  in  the  communications  constellation  s*  and  only  one  impulse  is 
allowed  in  each  symbol  period.  In  the  sparse  channel  identification  method  proposed 
by  [1],  estimation  of  the  ray  locations  is  performed  prior  to  estimation  of  ray  magni¬ 
tude.  The  ray  location  estimation  is  analogous  to  estimating  the  symbol  dither  and 
could  be  used  in  the  cooperative  receiver  MSE  approach  from  Section  3.2.21 
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IV.  Numerical  Simulation  Design 

In  this  chapter,  a  numerical  simulation  is  developed  to  investigate  the  LPI  sig¬ 
nal  design  from  Chapter  11111  The  model  includes  a  transmitter,  dispersive  noisy 
channels,  and  receivers  employing  several  equalization  techniques.  The  simulation 
calculates  the  experimental  bit  error  rate  for  a  transmitter,  channel,  receiver  set. 

The  Matlab®  model  is  described  in  Section  |4.1  and  the  associated  parameters 
are  outlined  in  Section  14.21  Modeling  limitations  and  assumptions  of  the  simulation 
are  discussed  in  Section  14.31  The  remainder  of  the  chapter  validates  the  components 
of  the  model  by  comparing  the  model  outputs  under  input  conditions  for  which  the 
result  is  predictable. 

4-1  Model  Overview 

The  end-to-end  communications  system  depicted  in  Figure  14. II  is  modeled  by 
a  discrete  numerical  Matlab®  simulation.  The  transmitter  block  converts  the  input 
bit  stream  into  a  baseband  waveform  representing  a  series  of  DPSK  symbols.  The 
waveform  is  a  discrete  sampling  of  the  analog  waveform  leaving  the  transmitter.  The 
channel  injects  white  Gaussian  noise  and  any  multipath  effects  into  the  system.  Fi¬ 
nally,  the  receiver  module  estimates  the  transmitted  symbols  from  the  received  signal 
(which  now  includes  channel  effects)  and  produces  the  output  bit  sequence. 

4-1.1  Transmitter.  The  transmitter  (Figure  14.21)  produces  the  sampled 
transmit  waveform  u(p)  from  the  input  bit  sequence.  The  bit  sequence  is  a  pseudo¬ 
random  binary  sequence  generated  by  the  Marsaglia  subtract  and  borrow  algorithm 
as  implemented  in  Matlab®  . 


{bin} 


{  bout } 


Figure  4.1:  Numerical  simulation  block  diagram 
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x(n)  Xup{p ) 


Figure  4.2:  Simulated  transmitter  block  diagram 


This  simulation  uses  Differential  Phase  Shift  Key  (DPSK)  modulation  to  encode 
the  input  bit  sequence  into  a  series  of  complex  symbols.  In  this  baseband  model, 
the  in-phase  and  quadrature  components  of  the  PSK  signal  are  represented  by  the 
real  and  imaginary  parts  of  the  symbol,  respectively.  In  the  complex  symbol  space, 
DPSK  constellations  are  circularly  symmetric  about  the  origin,  with  equal  magnitude 
symbols  and  constant  phase  difference  between  adjacent  symbols.  In  a  constellation 
with  M  symbols,  the  input  bit  stream  is  segmented  into  words  of  length  l  =  /og2(dT ) 
bits.  Each  /-bit  word  is  mapped  to  a  phase  offset  using  Gray  coding.  Per  (4.1)1 ,  the 
corresponding  symbol  is  found  by  adding  the  phase  offset  to  the  phase  of  the  previous 
symbol  in  the  series  [16:195]. 


I  si,  n=  1 

n)  =  < 

^  cos  {6n)  +  j  sin  (0n)  ,n>  1 

$n  =  &n- 1  + 


(4.1) 


Figure  14.31  shows  a  4-ary  DPSK  constellation  with  the  associated  phase  shift  mapping 
in  Table  14.11  In  this  example,  the  input  bit  stream  110001  results  in  the  symbol 
sequence 

[x(l),x(2),x(3),x(4)]  =  [si,  s3,  s3,  34]  ,  (4.2) 

where 


Si  =  cos  (0i)  +  j  sin  (0i) 
1  + 
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Figure  4.3:  4-ary  DPSK  Constellation 


Table  4.1:  4-ary  DPSK  Gray  code  mapping 


Word 

Phase  Shift,  A 6 

00 

0 

01 

7 r 

2 

11 

71 

10 

37 r 

2 
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Noise 

u(p) 


Channel 

Figure  4.4:  Simulated  channel  block  diagram 

The  Upsampling  Block  increases  the  sampling  frequency  by  inserting  an  integer 
number  of  zero-valued  samples  between  each  symbol.  After  upsampling,  the  symbols 
are  now  complex  impulses  spaced  P  samples  apart  (so  the  sampling  rate  is  P/Tsym 
samples  per  second).  Discrete  representations  of  continuous  time  signals  in  the  nu¬ 
merical  simulation  are  indexed  by  p. 

Symbol  timing  dither  to  create  ISI  is  implemented  by  the  Delay  Block.  The 
dither  is  realized  by  moving  each  impulse  some  integer  number  samples  according  to 
the  delay  code. 

Finally,  the  data  signal  consisting  of  complex  impulses  is  convolved  with  the 
pulse  shape  signal  f(p)  to  produce  the  transmitted  waveform. 

4-1.2  Channel.  The  simulated  physical  channel  is  a  combination  of  disper¬ 
sion  (h phys)  and  AWGN  [y{jp])  as  seen  in  Figure  4.41  Dispersion  (or  multipath)  effects 
are  included  by  convolving  the  transmitted  waveform  with  the  impulse  response  of  the 
physical  channel.  All  physical  channel  models  in  the  simulation  are  static.  Because 
this  is  a  baseband  simulation,  the  phase  distortion  that  could  be  caused  by  multipath 
carrier  signal  interference  is  not  directly  included  in  this  model.  To  indirectly  include 
phase  distortion,  the  coefficients  of  the  physical  channel  impulse  response  are  allowed 
to  be  complex. 
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Figure  4.5:  Simulated  receiver  block  diagram 


After  applying  the  multipath  effects,  AWGN  is  added  to  the  signal.  The  discrete 
noise  samples  are  complex,  with  both  the  real  and  imaginary  parts  of  each  noise 
samples  following  an  independent  Gaussian  distribution.  The  Marsaglia  ziggurat 
algorithm  in  Matlab®  is  used  to  produce  pseudorandom  Gaussian  samples.  The 
variance  of  both  the  real  and  imaginary  Gaussian  distributions  is  scaled  to  achieve  a 
desired  SNR  at  the  channel  output  (equivalent  to  the  input  of  the  receiver). 


4- 1. 3  Receiver.  The  simulated  receiver  model  performs  matched  filtering, 
downsampling,  post-processing,  and  demodulation  (Figure  14.5]).  The  matched  fil¬ 
ter  convolves  the  received  signal  with  the  same  pulse  shape  used  in  the  transmitter. 
Downsampling  of  the  filtered  signal  to  a  sampling  period  of  Tsym/ 2  results  in  two 
samples  per  symbol  period.  Although  the  simulation  is  expandable  to  downsample 
at  higher  rates,  fractionally  spaced  equalization  algorithms  are  restricted  to  two  sub¬ 
channels  in  this  thesis.  Receiver  models  using  symbol  spaced  data  use  only  the  first 
subchannel  y^(n).  Note  that  for  symbol  spaced  receiver  structures,  downsampling 
synchronization  significantly  impacts  the  SNR  at  the  downsampler  output.  This  sim¬ 
ulation  synchronizes  the  first  subchannel  to  the  point  where  the  impulse  response  of 
the  effective  channel  (including  the  effects  of  the  pulse  shape,  physical  channel,  and 
matched  filtering)  has  maximum  magnitude.  This  strategy  maximizes  the  SNR  for  the 
standard  transmitter  (no  pseudorandom  symbol  delay)  under  line-of-sight  conditions. 
When  the  pseudorandom  dither  transmitter  is  used,  the  first  subchannel  synchronizes 
to  the  average  location  of  the  maximum  effective  channel  impulse  response  within  the 
symbol  period. 
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Both  continuously  adaptive  and  periodically  adaptive  equalization  algorithms 
are  modeled  in  this  thesis.  To  directly  compare  the  bit  error  rates  of  continuously 
and  periodically  adaptive  methods,  the  equalizer  coefficients  are  held  constant  during 
data  collection.  The  coefficients  are  determined  by  training  the  equalizer  prior  to  each 
experiment  using  an  independent  series  of  bits. 

The  DPSK  demodulator  uses  the  phase  difference  between  the  current  complex 
estimated  symbol  and  the  previous  estimated  symbol  to  map  to  a  bit  string.  Because 
the  bit  estimation  is  based  only  on  the  phase  difference  between  successive  symbols, 
this  DPSK  demodulation  technique  is  not  affected  by  rotation  of  the  constellation 
caused  by  phase  ambiguity.  Sklar  calls  this  differentially  coherent  detection  of  differ¬ 
entially  encoded  PSK  modulation  [16:218].  For  binary  DPSK,  the  probability  of  bit 
error  is 

Pb  =  ,  (4-3) 

where  E b  is  the  information  signal  average  energy  per  bit  and  Na  is  the  noise  power 
spectral  density.  Eb  is  derived  from  the  average  energy  per  symbol  Es  and  the  number 
of  symbols  in  the  constellation  M\ 


Eb 


Es 

log2(M)' 


(4.4) 


Because  DPSK  demodulation  considers  both  the  current  symbol  and  the  previous 
symbol,  Es  includes  the  energy  of  both. 

For  differentially  coherent  detection  of  ilf-ary  DPSK,  the  approximate  proba¬ 
bility  of  symbol  error  is  [16:230] 


Pe(m)  “ 2Q  (7lrs inyii)  ■  (4-5) 

This  approximation  is  valid  only  for  large  Es/N0  ratios.  The  large  Es/N0  approxima¬ 
tion  leads  to  the  assumption  that  most  symbol  errors  result  in  misidentification  as  an 
adjacent  symbol,  producing  only  one  bit  error  per  symbol  error  due  to  Gray  coding. 
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In  this  manner,  the  probability  of  bit  error  for  M-ary  DPSK  is  estimated  by  [16:230] 


A(M)  « 


Pe(M) 

log 2(M)  ' 


(4.6) 


The  simulation  produces  an  output  bit  sequence  that  is  a  truncated  estimate  of 
the  input  bit  sequence.  This  truncation  is  due  to  the  convolution  operation  used  to 
apply  the  effective  channel  and  the  equalization  filter.  The  linear  convolution  output 
is  restricted  to  the  terms  for  which  both  signals  in  the  convolution  overlap  completely. 
In  this  manner,  the  convolution  of  signals  with  M  and  N  samples  produces  a  signal 
with  |  N  —  M |  +  1  samples.  The  net  effect  is  an  estimated  symbol  series  x(ri)  that  is 
shorter  than  the  transmitted  symbol  sequence  x(ri).  All  observed  bit  error  rates  in 
the  simulation  only  consider  the  bits  estimated  at  the  output  of  the  receiver,  ignoring 
the  leading  and  trailing  truncated  bits. 


4-2  Model  Parameters 

Many  transmitter,  channel,  and  receiver  parameters  are  variable  in  the  simula¬ 
tion.  The  large  number  of  parameters  creates  a  need  to  be  explicit  when  describing  a 
particular  simulation  result.  This  section  discusses  the  implementation  of  each  mod¬ 
ule  (transmitter,  channel,  and  receiver)  and  defines  the  relevant  parameters.  The 
simulation  parameters  are  summarized  in  Table  14121 

4-2.1  Global  Parameters.  In  numerical  modeling,  all  analog  signals  are  ap¬ 
proximated  by  discrete  signals.  The  fidelity  of  the  discrete  approximation  in  this 
simulation  is  determined  by  the  number  of  samples  per  symbol  period,  P.  All  exper¬ 
iments  presented  herein  used  64  samples  per  symbol  period  unless  otherwise  noted. 

4-2.2  Transmitter.  Within  the  Transmitter  module,  the  signal  constellation, 
pulse  shape,  and  symbol  delay  parameters  must  be  defined.  For  the  purposes  of  this 
simulation,  the  signal  constellation  and  pulse  shape  are  always  assumed  to  be  known 
by  the  receiver. 
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Table  4.2:  Simulation  parameters 


p 

Number  of  samples  per  symbol  period  used  to  model  ana¬ 
log  signals 

M 

Number  of  symbols  in  the  DPSK  constellation 

Pulse  Shape 

Pulse  shape  in  the  transmitter  and  the  corresponding 
matched  filter  in  the  receiver 

d 

Delay  code  sequence  used  to  determine  the  number  of 
samples  to  delay  each  symbol  in  the  transmitter 

Pmax 

Maximum  allowable  delay  in  the  sequence  d 

h 

LLphys 

Physical  channel  model 

SNR 

Signal-to-noise  ratio  as  measured  the  receiver  input 

Equalization 

Parameters 

Equalization  algorithm  and  associated  parameters 

0 

Number  of  samples  per  symbol  period  after  receiver 
downsampling  (Equivalent  to  the  number  of  subchannels 
in  the  multiple  output  model) 
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Table  4.3:  DPSK  Gray  code  mapping 


Binary 

4- ary 

8-ary 

Word 

Phase  Shift,  A 6 

Word 

Phase  Shift,  A 6 

Word 

Phase  Shift,  Ad 

0 

0 

00 

0 

000 

0 

1 

71 

01 

7 r 

2 

001 

7T 

4 

10 

37T 

2 

010 

37 r 

4 

11 

71 

Oil 

7 T 

2 

100 

7tt 

4 

101 

3n 

2 

110 

71 

111 

57 r 

4 

Binary,  4-ary,  and  8-ary  DPSK  Gray  coded  constellations  are  supported  by 
the  simulation  (Figure  14.6]) .  Each  constellation  is  circularly  symmetric  with  equal- 
magnitude  symbols,  satisfying  the  apriori  assumptions  of  all  the  equalization  tech¬ 
niques  discussed  in  Chapter  jl)  Table  14.31  gives  the  Gray  code  phase  shift  assignments 
for  each  DPSK  constellation. 

The  pulse  shapes  modeled  are  impulse,  square,  and  square  root  raised  co¬ 
sine  (SRRC),  plotted  in  Figure  14.71  The  impulse  pulse  shape  is  a  single  discrete 
impulse  with  unity  magnitude.  The  square  pulse  shape  is  a  rectangular  function  with 
width  equal  to  the  symbol  period  Tsym.  To  conserve  energy  through  the  Pulse  Shape 
block  in  the  transmitter,  the  amplitude  of  the  square  pulse  shape  is  1/Tsym. 

For  the  SRRC  pulse  shape,  the  roll-off  factor  (3  determines  the  “spread”  of  the 
shape  according  to  (]4.7j)  [6:225].  Because  the  SRRC  is  unbounded  in  time,  the  pulse 
must  be  truncated  in  any  physical  system.  When  the  SRRC  is  used  in  this  simulation, 
the  default  pulse  window  is  ten  symbol  periods  wide,  centered  at  t  —  0  in  the  SRRC 


4-9 


Figure  4.6:  (a)  Binary  DPSK  Constellation  (b)  4-ary  DPSK  Constellation  (c)  8-ary 

DPSK  Constellation 
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pulse  shape  equation: 


/(*) 


1  sm[n(l—f3)t/T8ym]  +  (4:(3t/TSym)  cos[7r(l +/3)t/Tsym\ 

yjTsym  (jf"t/Tsym)  [l  —  (^f3t/TSym)  ] 


[1  ~P  +  (4/3/vr)] , 


yyTsyrn 

-  (1  +  f)sin(ih)  +  (1-f)cos(§ 


2T„, 


t  =  0 


(4.7) 


To  implement  the  symbol  timing  dither,  the  Delay  Block  delays  each  symbol 
impulse  up  to  P  —  1  samples.  The  delay  code  d  is  a  pseudorandom  series  of  integers 
that  defines  the  particular  delay  of  each  symbol  (4.8).  Being  a  finite  pseudorandom 
series,  the  pattern  of  delays  repeats  every  N(i  symbols,  where  Nd  is  the  length  of  the 
delay  code.  The  maximum  allowable  delay  must  be  set  to  less  than  P  —  1  samples 
to  prevent  adjacent  symbols  from  being  transmitted  in  the  wrong  order  (]4.9j).  The 
distribution  of  integer  delays  in  the  random  delay  sequence  is  always  uniform  in  this 
simulation.  To  model  a  standard  transmitter,  the  maximum  allowable  delay  is  simply 
restricted  to  zero  samples. 


d  =  (d(l),d(2),...,d(iVd)} 
d{n)  G  I[0,  Pmax\i  0  ^  Pmax  ^  P  1 


(4.8) 

(4.9) 


4-2.3  Channel.  The  simulated  physical  channel  is  a  combination  of  disper¬ 
sion  and  AWGN.  The  dispersion  is  applied  via  convolution  with  the  physical  channel 
impulse  response  per  Section  14.11  Although  any  discretized  impulse  response  can  be 
inserted  into  the  simulation,  this  thesis  characterizes  receiver  performance  using  a 
line-of-sight,  lossless  channel  with  a  delay  of  ten  samples. 

AWGN  is  added  to  the  signal  after  applying  the  physical  channel  impulse  re¬ 
sponse.  The  variance  of  the  noise  samples  can  be  specified  either  by  SNR  or  Eb/N0  at 
the  receiver  input.  The  actual  simulated  SNR  is  then  estimated  by  taking  the  ratio 
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Figure  4.7:  (a)  SRRC  pulse  shape,  /3  —  0  (b)  SRRC  pulse  shape,  (3  =  0.25  (c)  Square 

pulse  shape  (d)  Impulse  pulse  shape 
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of  mean  squares  of  the  information  signal  z(p)  and  noise  signal  v(p)  according  to 

SNR  «  .  (4.10) 

EHp)\2 

In  the  line-of-sight  case  with  no  multipath  effects  and  no  symbol  dither,  the  inter¬ 
symbol  interference  is  negligible  (at  least  for  the  pulse  shapes  used  in  this  thesis),  so 
the  signal  power  in  the  SNR  calculation  can  be  estimated  directly  from  the  simula¬ 
tion  just  prior  to  adding  the  noise  in  Figure  14.41  When  ISI  is  present  in  the  signal, 
some  clarification  is  necessary.  For  the  purposes  of  this  document,  the  signal  power  is 
still  measured  directly  from  the  simulated  signal  prior  to  adding  noise  (after  applying 
multipath  effects)  and  noise  power  is  still  measured  from  the  AWGN  signal  v(p).  In 
practice,  however,  the  ISI  due  to  multipath,  symbol  dither,  and/or  non-ideal  pulse 
shapes  can  be  treated  as  an  additional  noise  source.  The  contribution  of  ISI  to  bit 
errors  is  dependent  on  the  quality  of  equalization,  with  perfect  equalization  eliminat¬ 
ing  ISI  completely.  Hence,  all  SNR  calculations  in  this  thesis  do  not  consider  ISI  as 
a  source  of  noise  power  unless  explicitly  stated. 

4.2.4  Receiver.  The  receiver  parameters  are  different  for  each  receiver  struc¬ 
ture  modeled  in  the  simulation.  Table  14.41  lists  the  algorithms  and  parameters  asso¬ 
ciated  with  each.  The  optimal  set  of  parameters  for  each  receiver  model  is  highly 
dependent  on  the  transmitter  and  channel  parameters.  Optimization  of  a  receiver 
parameter  set  is  not  trivial  due  to  the  number  of  permutations  of  parameter  combi¬ 
nations.  In  lieu  of  exhausting  all  possible  receiver  configurations,  the  parameters  for 
each  experiment  are  chosen  by  making  reasonable  assumptions  about  the  receivers 
knowledge  of  the  transmitter  and  the  channel. 

For  algorithms  using  statistical  averaging,  N  is  the  number  of  symbols  output 
by  the  transmitter  for  the  purpose  of  estimating  the  equalization  hlter(s).  The  actual 
number  of  symbols  available  to  the  algorithm  in  the  receiver  is  somewhat  less  than 
N  due  to  truncation  of  the  leading  and  trailing  samples.  This  truncation  omits  the 
received  samples  that  reflect  initiation  and  termination  transients.  For  most  results 
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Table  4.4:  Equalization  algorithm  parameters 


Basic 

N/A 

No  parameters  necessary 

FS  CMA 

h 

Gradient  descent  step  size 

M 

Number  of  taps  in  the  equalizer  (equalizer  length) 

N 

Number  of  symbols  used  to  train  the  equalizer 

Subspace 

L 

Estimate  of  the  length  of  the  effective  channel  (in  symbol 
periods) 

M 

Number  of  taps  in  the  equalizer  (equalizer  length) 

N 

Number  of  symbols  used  to  train  the  equalizer 

FS  LMS  Bank 

/i 

Gradient  descent  step  size 

M 

Number  of  taps  in  the  equalizer  (equalizer  length) 

5 

Impulse  response  delay  to  which  the  equalizer  converges 
(in  symbol  periods) 

N 

Number  of  symbols  available  for  equalizer  estimation 

Modified  MSE 

6 

Impulse  response  delay  to  which  the  equalizer  converges 
(in  symbol  periods) 

L 

Estimate  of  the  length  of  the  effective  channel  (in  symbol 
periods) 
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presented  in  this  thesis,  N  —  5,  000  symbols  which  is  be  large  enough  for  algorithms 
based  on  higher-order  statistics  to  reach  steady-state.  This  is  required  to  reflect  the 
“best  case”  bit  error  performance.  It  should  be  acknowledged  that  this  approach 
renders  the  superior  convergence  speed  of  the  subspace  decomposition  algorithm  a 
moot  point  in  this  simulation. 

The  gradient  descent  algorithms,  fractionally  spaced  CMA  and  the  bank  of 
fractionally  spaced  LMS  equalizers,  require  an  iteration  step  size  /i.  The  magnitude 
of  the  step  size  must  strike  a  balance  between  maximizing  convergence  speed  and 
minimizing  steady-state  error  (see  Section  |2.4.2j).  For  this  simulation,  /i  has  been 
tuned  for  each  set  of  transmitter /channel  parameters  using  emperical  simulation. 

The  number  of  taps  in  the  equalization  filter,  M,  must  be  defined  for  the  trained 
algorithms.  For  the  FS  CMA  and  FS  LMS  bank  receivers,  the  number  of  taps  is 
chosen  to  satisfy  the  minimum  length  condition  in  (12.301)  for  (potentially)  perfect 
equalization,  (see  Section  12.4.3).  The  maximum  effective  channel  length  is  L  =  20 
symbol  periods,  so  M  =  40  unless  otherwise  noted  for  these  two  receivers. 

The  subspace  decomposition  algorithm  equalization  solution  is  unstable  when 
the  effective  channel  length  estimate  L  and  equalizer  length  M  are  overestimated.  To 
prevent  the  potential  noise  amplification  that  occurs  when  this  happens,  these  two 
parameters  are  tuned  empirically  by  observing  the  output  bit  error  rate  for  a  range 
of  parameter  sets  for  each  transmitter /channel  configuration. 

For  the  FS  LMS  bank  and  modified  MSE  algorithms,  the  6  parameter  establishes 
the  difference  between  the  indexing  of  the  input  bit  stream  and  the  estimated  output 
bit  stream.  For  an  unknown  channel,  <5  is  often  determined  by  minimizing  bit  error 
rate  over  a  range  of  5  values.  In  simulation,  pulse  shape  duration,  physical  channel 
impulse  response,  and  the  matched  filter  are  known  and  combined  to  calculate  the 
effective  channel  impulse  response  length  L.  With  this  information,  each  matched 
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filter  sample  y^(n)  is  known  to  be  a  function  of  the  previous  L  symbols  according  to 


y^l\n)  —  F  [x(n  —  L  +  1),  x{n  —  L  +  2), x(n)\  .  (4-11) 


In  a  similar  manner,  the  estimated  symbol  x(ri)  is  a  function  of  previous  M/O  matched 
filter  samples  in  each  subchannel  y^\n).  Combining  the  effective  channel  impulse 
response  length  and  number  of  equalization  filter  taps  allows  an  estimated  symbol  to 
be  a  function  of  a  specific  set  of  transmitted  symbols: 


x(n)  =  F 


x  n 


M 

~o 


L  +  2 


,  x  \  n 


M 

~o 


L  +  l 


,x(n 


(4.12) 


Combining  (]4.11j)  and  (14.1211  provides  the  natural  bounds  on  6: 


&opt  £  s  0, 1, L  + 


M 

~o 


(4.13) 


This  thesis  studies  line-of-sight  channels  with  a  known  transmission  delay.  For 
an  ideal  linc-of-sight  channel,  the  effective  channel  impulse  response  is  the  convolution 
of  the  pulse  shape  and  the  matched  filter  with  some  delay,  which  is  symmetric  about 
the  impulse  response  maximum.  In  the  particular  linc-of-sight  channel  used  in  this 
thesis,  the  transmission  delay  is  ten  samples,  or  17%  of  a  symbol  period,  so  a  good 
estimate  for  5opt  is 


(4.14) 


Unless  otherwise  stated,  this  is  the  value  of  5  used  in  this  thesis. 


For  the  modified  MSE  receiver,  the  number  of  equalizer  filter  taps  M  is  a  func¬ 
tion  of  the  estimated  channel  length  L  according  to 


M  —  OL 


(4.15) 


For  the  line-of-sight  channel  used  in  this  thesis,  the  actual  channel  length  L  is  used 
as  the  estimated  channel  length  parameter  for  the  modified  MSE  receiver. 
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4-3  Assumptions 

The  simulation  described  in  this  chapter  makes  several  assumptions  and  simpli¬ 
fications  in  the  communications  system  model.  These  limitations  are  detailed  below. 

By  using  a  baseband  simulation,  the  effects  of  carrier  waveform  phase  interfer¬ 
ence  between  multipath  reflections  is  not  captured  explicitly  in  the  multipath  model. 
Instead,  the  potential  phase  interference  effects  of  the  carrier  waveform  are  accounted 
for  by  the  baseband  multipath  channel  coefficients  in  the  simulation.  These  coeffi¬ 
cients  are  allowed  to  be  complex  to  model  the  impact  of  the  reflection  on  the  phase 
of  the  carrier  signal. 

Doppler  effects  are  not  modeled  in  this  baseband  simulation.  The  frequency  of 
the  notional  carrier  waveform  is  assumed  to  be  equal  for  each  multipath  reflection  at 
the  receiver  input. 

The  receiver  knows  the  time  when  the  impulse  response  of  the  combined  chan¬ 
nel/pulse  shape  is  at  a  maximum  and  synchronizes  the  primary  channel  to  sample  at 
this  point.  This  synchronization  is  critical  for  symbol  spaced  receivers.  Conceivably,  a 
blind  detector  could  estimate  this  sampling  synchronization  via  energy  detection  and 
averaging  over  many  symbols.  When  the  transmitter  applies  pseudorandom  delays 
to  each  symbol,  the  point  within  each  symbol  period  for  which  the  effective  channel 
impulse  response  is  at  a  maximum  in  the  receiver  is  no  longer  the  same  for  each  sym¬ 
bol.  In  this  case,  the  receiver  synchronizes  to  the  average  time  at  which  the  impulse 
response  is  at  a  maximum.  Again,  this  synchronization  is  deterministic  within  the 
simulation,  so  the  receiver  synchronization  is  assumed  to  be  perfect  and  stable. 

The  modified  MSE  receiver  requires  knowledge  of  the  additive  channel  noise 
power.  This  simulation  uses  sample  mean  and  variance  of  the  noise  signal  to  directly 
determine  the  noise  power.  In  reality,  the  noise  signal  is  not  available  independent 
from  the  information  signal  at  the  receiver  input,  so  another  noise  estimation  tech¬ 
nique  would  be  required.  This  could  be  as  simple  as  measuring  the  power  at  the 
receiver  input  when  the  transmitter  is  not  present. 
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Figure  4.8:  Simulated  bit  error  rates  for  DPSK  in  a  LOS  channel  compared  to  the 
analytical  probability  of  bit  error. 

4-4  Simulation  Verification 

This  section  presents  incremental  simulation  results  to  validate  the  simulation 
components.  To  demonstrate  the  modulator  and  demodulator  models,  the  bit  error 
rates  produced  via  simulation  under  line-of-sight  conditions  are  compared  to  theo¬ 
retical  DPSK  performance.  The  effects  of  multipath  models  are  demonstrated  by 
analyzing  the  impact  on  bit  error  rate  and  examining  the  received  symbols  in  the 
signal  space.  The  equalizers  are  then  inserted  into  the  receiver  model  to  demonstrate 
their  multipath  mitigation  capabilities. 

4-4-1  DPSK  Simulation  Validation.  Figure  14.8  compares  the  analytical 
probability  of  bit  error  for  DPSK  constellations  in  the  presence  of  AWGN  to  the 
bit  error  rates  observed  in  simulation  using  an  ideal  line-of-sight  physical  channel. 
Over  the  Eb/N0  range  considered  here,  the  difference  between  analytical  bit  error 
probability  and  the  simulated  result  is  not  statistically  significant.  All  error  bars 
on  plots  in  this  thesis  represent  95%  confidence  intervals  produced  by  Monte  Carlo 
simulation  [10:266].  An  ideal  square  pulse  shape  is  used  in  the  transmitter  to  prevent 

ISI. 
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Figure  4.9:  Simulated  bit  error  rates  for  binary  DPSK  in  a  LOS  channel.  The 

truncated  SRRC  shape  introduces  ISI  that  degrades  the  receiver  performance. 

In  practice,  the  SRRC  pulse  shape  is  more  practical  than  an  ideal  square  wave 
because  of  bandwidth  considerations.  An  infinite  length  SRRC  pulse  shape  and 
matched  filter  pair  results  in  no  ISI,  like  the  square  wave  pulse  shape.  When  the 
SRRC  pulse  shape  is  (necessarily)  truncated,  typically  to  a  total  width  of  10 Tsym  in 
this  simulation,  the  symbol-spaced  samples  at  the  matched  filter  output  are  subject 
to  ISI,  degrading  the  bit  error  rate.  Figure  [4AH  demonstrates  the  effect  of  ISI  due  to 
SRRC  windowing.  Note  that  for  equal  size  pulse  shape  windows,  the  bit  error  rate 
degradation  is  most  severe  for  small  (3  (wide  SRRC  shape).  Figure  14.101  shows  that 
small  (3  results  in  large  sidelobes  in  the  effective  channel  impulse  response  with  zero 
crossings  that  do  not  synchronize  to  the  rate  1/P  downsampling. 

4-4 Multipath  Interference  Validation.  A  simple  sparse  physical  channel 
model  shown  in  Figure  14.11(a)  is  used  to  demonstrate  multipath  modeling  in  the 
simulation.  The  effect  of  the  channel  on  the  location  of  received  symbols  in  the  signal 
space  of  a  QPSK  modulated  waveform  (at  Eb/No  =  9  dB  measured  at  the  matched 
Liter  input)  is  given  in  Figure  [4.11(b) .  An  SRRC,  (3  =  0.25  pulse  shape  and  matched 
Liter  are  in  use  here.  After  differentially  coherent  demodulation,  the  observed  BER  is 
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Effective  Channel,  SRRC  p  =  0,  LOS  Effective  Channel,  SRRC  p  =  0.25,  LOS 


Figure  4.10:  Effective  channel  impulse  responses  (pulse  shape  and  matched  filter). 
The  receiver  synchronizes  the  first  symbol  spaced  subchannel  to  the  impulse  response 
maximum,  (a)  (3  =  0.  (b)  (3  —  0.25 

0.24  vice  the  BER  of  0.0088  observed  in  simulation  for  the  line-of-sight  channel  at  this 
bit  energy  to  noise  ratio.  Note  that  the  imaginary  component  of  the  physical  channel 
impulse  response  produces  a  rotation  of  the  constellation  in  the  symbol  space.  The 
BER  degradation  due  to  this  multipath  channel  is  summarized  over  a  range  of  E^/N0 
ratios  in  Figure  IT.  121  Over  this  range,  the  severe  ISI  caused  by  the  multipath  channel 
dominates  the  accuracy  of  the  demodulator  bit  estimation;  the  8  dB  gain  in  Eb/N0 
only  improves  the  BER  from  0.28  to  0.23. 

4-4-3  Equalizer  Validation.  To  demonstrate  the  ISI  mitigation  capability 
of  fractionally  spaced  CMA,  fractionally  spaced  LMS,  and  subspace  decomposition, 
each  equalizer  is  applied  transmitter  and  channel  described  in  the  previous  section. 
Each  equalizer  uses  N  —  5,  000  transmitted  symbols  for  training  before  the  equalizer 
output  symbols  are  demodulated  for  BER  calculations.  After  demodulating  5,000 
symbols,  the  equalizers  are  reset  and  trained  again,  analogous  to  periodically  adaptive 
equalization.  For  fractionally  spaced  CMA  and  LMS,  the  step  size  /i  is  selected 
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Physical  Channel  Impulse  Response 


Received  Symbols,  No  Equalization 


Figure  4.11:  (a)  Impulse  response  of  a  sparse  multipath  channel,  (b)  Received 

samples  in  the  signal  space  for  QPSK  modulation  using  an  SRRC  (3  =  0.25  pulse 
shape  at  Eb/Na  =  9  db. 


Figure  4.12:  ISI  due  to  the  sparse  multipath  channel  model  dominates  the  observed 
BER, 
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Figure  4.13:  Applying  equalizer  filters  at  the  matched  filter  output  improves  BER 
performance  for  the  multipath  channel.  The  4-ary  PSK  communications  system  is 
using  an  SRRC,  (3  =  0.25  pulse  shape  and  matched  filtering. 

experimentally  to  produce  near-optimal  convergence  over  the  Eb/N0  range  of  interest 
in  the  sparse  channel  environment  with  the  5,000  symbols  transmitted  for  training. 
The  effective  channel  length  L  and  the  number  of  equalizer  filter  taps  M  are  also 
determined  empirically. 

Figure  [4.131  compares  the  bit  accuracies  achieved  by  CMA,  LMS,  and  subspace 
decomposition  equalization  algorithms  to  the  BER  observed  without  equalization. 
CMA  and  LMS  error  rates  behave  similarly  across  the  Eb/N0  range  shown  here.  None 
of  the  equalizer  outputs  approach  the  analytical  probability  of  bit  error  achievable  in 
the  absence  of  multipath  interference.  The  subspace  method  shows  only  a  slight 
improvement  over  the  unequalized  case  for  this  multipath  channel. 
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V.  Simulation  Results 


This  chapter  presents  numerical  simulation  results  demonstrating  the  LPI  charac¬ 
teristics  of  the  pseudorandom  dither  communications  signal  described  in  Chap¬ 
ter  HU  The  performance  of  both  cooperative  receivers  and  non-cooperative  receivers 
are  quantified  by  the  observed  bit  error  rate  after  signal  demodulation.  Section  15711 
characterizes  the  behavior  of  the  cooperative  receivers  (modified  MSE  and  LMS  bank). 
The  non-cooperative  receiver  simulation  results  are  presented  in  Section  15.2.1  The  co¬ 
operative  and  non-cooperative  receivers  are  compared  in  Section  15151  in  an  attempt  to 
select  a  set  of  operating  parameters  that  yields  LPI  characteristics  with  little  penalty 
to  the  cooperative  receiver  link.  Finally,  the  bandwidth  impact  of  the  dither  strategy 
on  DPSK  modulation  is  considered  in  Section  15.41 

5.1  Cooperative  Receiver  Performance 

Results  for  the  modified  MSE  deconvolution  and  adaptive  LMS  equalizer  bank 
receiver  structures  are  presented  below.  The  performance  of  each  cooperative  re¬ 
ceiver  structure  is  demonstrated  using  the  same  set  of  analyses.  First,  the  bit  error 
rates  (BER)  oberved  for  each  cooperative  receiver  using  traditional  binary,  4-ary,  and 
8-ary  DPSK  signalling  schemes  are  presented  (using  a  line-of-sight  (LOS),  AWGN 
channel).  This  result  is  compared  to  the  analytical  probability  of  bit  error  (Pb)  for 
each  particular  signalling  scheme  to  identify  any  potential  BER  penalty  associated 
with  the  receiver  structure  for  standard  DPSK  modulated  signals.  Note  that  standard 
M-ary  DPSK  modulation  is  equivalent  to  using  the  PDW  with  the  symbol  dither  set 
to  zero  seconds  for  each  symbol.  Next,  the  BER  behavior  of  each  receiver  is  exam¬ 
ined  with  respect  to  the  pseudorandom  dither  waveform  (PDW)  dither  range.  BER 
results  for  varying  SRRC  pulse  shape  roll-off  factor  (/3)  and  symbol  dither  range  are 
contrasted  with  the  goal  of  selecting  a  set  of  waveform  parameters  for  each  coop¬ 
erative  receiver  structure  that  allows  for  accurate  detection  of  the  PDW  signal  and 
significant  detection  degradation  in  non-cooperative  receivers.  This  set  of  parameters 
for  the  PDW  is  used  to  study  BER  behavior  for  a  range  of  Eb/Na. 
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Figure  5.1:  The  modified  MSE  cooperative  receiver  accuracy  is  statistically  similar 
to  the  analytical  probability  of  bit  error. 

5.1.1  Modified  MSE  Deconvolution.  For  a  LOS  channel,  the  observed  BER 
using  the  modified  MSE  deconvolution  receiver  structure  matches  the  analytical  bit 
error  probability  predicted  by  (14.31).  (14.5).  and  (4.61)  for  input  Eb/Na  ratios  shown 
in  Figure  15711  The  error  bars  on  the  observed  BER  data  points  in  this  plot  and  all 
others  in  this  chapter  represent  95%  confidence  intervals  produced  by  Monte  Carlo 
simulation  [10:266].  In  addition,  the  complexity  of  the  modified  MSE  receiver  is 
significantly  reduced  when  the  symbol  dither  is  restricted  to  zero  because  only  a 
single  convolution  matrix  is  required,  and  the  inverse  need  only  be  calculated  a  single 
time.  Thus,  equalization  using  the  modified  MSE  receiver  in  a  LOS,  AWGN  channel 
is  essentially  equivalent  to  a  simple  matched-filter  receiver  without  post-processing. 

Figure  15121  shows  the  modified  MSE  receiver  error  rate  behavior  for  varying  sym¬ 
bol  dither  ranges.  The  results  reflect  a  simulated  LOS,  AWGN  channel  and  Eb/N0  =  9 
dB  at  the  receiver  input.  As  the  symbol  dither  range  increases  to  a  full  symbol  pe¬ 
riod,  the  severity  of  the  ISI  increases,  causing  the  observed  BER  to  degrade.  The 
dither  threshold  at  which  the  BER  diverges  from  the  zero  dither  level  is  dependent 
on  the  pulse  shape.  Specifically,  the  more  narrow  the  pulse  shape  (corresponding  to 
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Figure  5.2:  Modified  MSE  receiver  performance  as  a  function  of  the  dither  range 

using  a  LOS  channel  at  Eb/N0  =  9  db.  (a)  Binary  DPSK  modulation  (b)  4-ary  DPSK 
modulation 
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Figure  5.3:  Applying  the  bank  of  LMS  equalizers  in  a  LOS  environment  increases 
the  BER  with  respect  to  the  demodulation  of  the  unequalized  matched  filter  output 
due  to  misadjustment. 

larger  j3),  the  larger  the  dither  that  can  be  tolerated  before  demodulation  accuracy 
is  degraded  significantly.  This  directly  follows  from  the  lower  1S1  contribution  due  to 
neighboring  symbols  for  narrow  pulse  shapes. 

5.1.2  Adaptive  LMS  Equalizer  Bank.  Unlike  the  modified  MSE  receiver, 
the  LMS  bank  of  equalizers  receiver  structure  produces  a  BER  slightly  worse  than 
the  analytical  probability  of  bit  error  for  standard  DPSK  modulation  schemes,  as 
seen  in  Figure  15.31  As  discussed  in  Section  |2.4.2[  the  step  size  ji  and  noise  prevent 
the  LMS  equalizers  in  the  receiver  from  converging  to  an  exact  solution.  This  effect 
can  be  mitigated  by  adaptively  adjusting  the  LMS  filter  lengths  via  channel  length 
estimation  and  decreasing  the  step  size  as  the  filter  converges,  which  are  beyond  the 
scope  of  this  thesis. 

The  correlation  between  symbol  dither  range  and  observed  BER  is  shown  in 
Figure  15.41  The  distinction  between  the  SRRC  f3  parameters  is  less  pronounced 
than  for  the  modified  MSE  receiver  because  misadjustment  dominates  the  minor  1SI 
variation  due  to  SRRC  (3.  As  with  the  modified  MSE  receiver,  1S1  due  to  symbol 
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Figure  5.4:  LMS  bank  receiver  performance  as  a  function  of  the  dither  range  using 

a  LOS  channel  at  Eb/N0  =  9  db.  (a)  Binary  DPSK  modulation  (b)  4-ary  DPSK 
modulation 
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dither  dominates  the  BER  result  as  the  dither  range  increases  to  an  entire  symbol 
period. 

5.2  Non-Cooperative  Receiver  Performance 

This  section  characterizes  the  performance  of  the  non-cooperative  receivers  de¬ 
scribed  in  Section  3.31  The  observed  error  rates  of  each  receiver  are  analyzed  in  a 
simulated  LOS,  AWGN  channel. 

5.2.1  Basic  Receiver.  The  basic  receiver  does  not  perform  any  ISI  mitiga¬ 
tion,  simply  demodulating  the  samples  of  the  matched  filter  output.  This  architecture 
produces  the  observed  BER  in  Figure  [5751  For  both  binary  and  4-ary  DPSK  modula¬ 
tion,  the  BER  diverges  from  the  the  analytical  P&  at  a  uniform  dither  range  of  0.2 Tsym. 
The  BER  appears  to  be  approaching  an  upper  limit  at  as  dither  range  increases.  Also, 
the  impact  of  SRRC  pulse  shape  truncation  is  more  pronounced  in  this  receiver  struc¬ 
ture  than  in  either  cooperative  receiver,  indicating  that  the  two  cooperative  receivers 
are  mitigating  this  ISI  source  as  well  as  the  symbol  dither  ISI. 

5.2.2  Blind  Equalizer  Receiver.  The  observed  BER  for  the  receiver  em¬ 
ploying  FS-CMA  equalization  in  a  LOS  environment,  shown  in  Figure  15.61  is  not 
significantly  improved  relative  to  the  basic  receiver  error  rates  in  Figure  [5751  In  fact, 
for  many  parameter  combinations,  FS-CMA  equalization  increases  the  bit  error  rate 
due  to  misadjustment  in  the  FS-CMA  steady-state.  It  should  also  be  noted  that  the 
particular  sampling  synchronization  technique  employed  in  this  simulation  minimizes 
the  benefit  of  fractionally  spaced  equalization  techniques  over  the  symbol  spaced  sam¬ 
pling  in  the  basic  receiver.  Interestingly,  the  ISI  due  to  SRRC  truncation,  which  is 
most  severe  for  small  /3,  is  not  distinctive  in  the  error  rates  for  low  dither  ranges, 
another  symptom  of  misadjustment.  The  same  phenomenon  can  be  seen  to  a  lesser 
degree  in  the  LMS  equalizer  bank  result.  In  both  cases,  the  (3  —  0  error  curve  is  lower 
than  that  of  the  basic  receiver,  indicating  that  some  equalization  of  the  ISI  due  to 
pulse  shape  truncation  is  occurring  when  ISI  due  to  symbol  dither  is  minimal. 
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Figure  5.5:  Basic  receiver  performance  as  a  function  of  the  dither  range  using 

a  LOS  channel  at  Eb/N0  =  9  db.  (a)  Binary  DPSK  modulation  (b)  4-ary  DPSK 
modulation 


5-7 


-1 


10  r 


-2 


DC 

LU 

CD 


10  - 


10 


-4 


P  =  0 
p  =  0.25 
(3  =  0.5 
Analytical  P 


io-3h . 


r 

Jr" 


-r 


0.1 


0.2 


0.3 


0.4  0.5 


0.6 


0.7 


0.8 


0.9 


Dither  range  as  a  fraction  of  T 


sym 


(a) 


(b) 


Figure  5.6:  FS-CMA  equalizer  receiver  performance  as  a  function  of  the  dither 

range  using  a  LOS  channel  at  Eb/N0  =  9  db.  (a)  Binary  DPSK  modulation  (b)  4-ary 
DPSK  modulation 
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The  error  rates  observed  in  simulation  of  the  subspace  equalization  receiver 
(Figure  15.7])  indicate  that  subspace  equalization  is  not  suited  for  mitigation  of  IS1 
due  to  symbol  dither  in  a  LOS  channel.  All  of  the  BER  curves  are  higher  than  their 
equivalents  in  the  basic  receiver.  Even  when  the  dither  is  removed,  the  inversion  of 
the  correlation  matrix  within  the  subspace  channel  identification  algorithm  is  poorly 
conditioned,  producing  the  observations  in  Figure  [573 

5.2.3  Symbol  Dither  Estimation  Receiver.  The  symbol  location  estimation 
algorithm  described  in  Section  3.3.3  produces  poor  symbol  identification  estimates, 
as  seen  in  Figure  15.81  All  bit  error  rates  shown  here  fall  between  7.9%  and  29%. 
Differentiation  between  SRRC  (3  values  is  only  apparent  when  the  symbol  dither  is 
allowed  to  vary  over  0.5 Tsym  or  greater.  The  improvement  in  BER  for  increasing 
symbol  dither  is  contrary  to  intuition  and  the  behavior  of  all  other  receivers.  This 
anomaly  is  most  likely  an  artifact  of  the  symbol  boundary  definition  within  the  receiver 
simulation  which  biases  the  error  in  the  estimated  symbol  locations  produced  by  the 
algorithm  for  small  dither  ranges.  Clearly,  a  more  rigorous  estimation  algorithm  is 
required  to  make  this  receiver  structure  viable. 


5-9 


-1 


10  r 


-2 


DC 

LU 

CD 


10  - 


10-3, 


10 


-4 


-P  =  0 
-(3  =  0.25 
-(3  =  0.5 
—  Analytical  P 


ie  : 


i: 


:: 


... 


-a: : 


0.1 


0.2 


0.3 


0.4  0.5 


0.6 


0.7 


0.8 


0.9 


Dither  range  as  a  fraction  of  T 


sym 


(a) 


(b) 


Figure  5.7:  Subspace  equalizer  receiver  performance  as  a  function  of  the  dither 

range  using  a  LOS  channel  at  Eb/N0  =  9  db.  (a)  Binary  DPSK  modulation  (b)  4-ary 
DPSK  modulation 
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Figure  5.8:  Dither  estimation  receiver  performance  as  a  function  of  the  dither  range 
for  binary  DPSK  modulation  using  a  LOS  channel  at  Ef,/N0  =  9  db. 


5.3  LPI  Signal  Parameter  Selection 

This  section  explores  the  potential  suitability  of  applying  the  pseudorandom 
dither  strategy  with  DPSK  modulation  as  an  LPI  signalling  scheme.  The  goal  is  to 
find  a  transmitter/cooperative  receiver  pair  that  1)  yields  a  BER  near  optimal  DPSK 
performance  and  2)  results  in  poor  BER  for  the  non-cooperative  receivers.  Several 
promising  parameter  sets  are  selected  using  the  cooperative  and  non-cooperative  per¬ 
formance  curves  presented  earlier  in  this  chapter.  The  relative  observed  BER  among 
the  cooperative  and  non-cooperative  receivers  suggests  that  symbol  dither  can  success¬ 
fully  be  employed  to  impede  exploitation  of  the  signal  information  by  non-cooperative 
receivers  while  preserving  the  cooperative  communication  link. 

To  select  a  set  of  signal  parameters,  the  performance  curves  are  grouped  by 
modulation  type  and  pulse  shape  in  Figure  f5T91  through  Figure  [5.141  In  these  figures, 
cooperative  receivers  BER  are  denoted  by  solid  lines  while  non-cooperative  receiver 
results  are  represented  by  dashed  lines.  The  greatest  bit  error  difference  (in  a  loga¬ 
rithmic  sense)  between  a  cooperative  receiver  and  the  most  accurate  non-cooperative 
receiver  occurs  for  the  binary  DPSK  modulated  signal  using  an  SRRC  pulse  shape 
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Figure  5.9:  Receiver  performance  comparison  for  binary  DPSK,  SRRC  (3  =  0  at 
Eb/N0  =  9  db. 


with  (3  =  0.25  and  a  dither  distribution  range  of  0.6 Tsym.  At  this  operating  point, 
found  in  Figure  15.101  the  modified  MSE  cooperative  receiver  yields  BER  =  5.4  •  10~4 
versus  the  best  non-cooperative  receiver  at  BER  =  8.9  •  10~3.  Bit  errors  occur  16.4 
times  as  frequently  in  the  non-cooperative  receiver  than  in  the  cooperative  receiver 
at  Eb/N0  =  9  dB. 

The  BER  improvement  between  cooperative  and  non-cooperative  receivers  is  not 
as  high  for  4-ary  DPSK  modulation  at  Eb/N0  =  9  dB,  largely  because  the  analytical 
4-ary  DPSK  BER  at  Eb/N0  =  9  dB  is  higher  than  the  binary  DPSK  BER.  Thus, 
the  observed  BER  for  the  different  receivers  is  dominated  by  AWGN  vice  receiver 
efficiency.  For  4-ary  DPSK  modulation  with  SRRC  pulse  shape  (3  =  0.25  and  symbol 
dither  range  of  0.6 T3ym,  the  non-cooperative  BER  is  only  7.4  times  worse  than  the 
modified  MSE  cooperative  receiver  BER,  as  seen  in  Figure  [5.131 

Fixing  the  symbol  dither  range  at  0.6Tsym  and  varying  Eb/N0  for  binary  and 
4-ary  constellations  produces  the  performance  results  shown  in  Figure  [5.151  and  Fig¬ 
ure  15.161  Between  the  two  constellations,  the  4-ary  scheme  produces  greater  separa¬ 
tion  between  cooperative  and  non-cooperative  receiver  performance.  For  example,  to 
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Figure  5.10:  Receiver  performance  comparison  for  binary  DPSK,  SRRC  {3  =  0.25 
at  Eb/N0  =  9  db. 


Figure  5.11:  Receiver  performance  comparison  for  binary  DPSK,  SRRC  /3  =  0.50 
at  Eb/N0  =  9  db. 
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Eb/N0  =  9  db. 


Eb/N0  =  9  db. 
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Figure  5.14:  Receiver  performance  comparison  for  4-ary  DPSK,  SRRC  (3  =  0.50  at 
Eb/N0  =  9  db. 

operate  at  BER  =  1.0  •  10~2  non-cooperative  receivers  require  ~  3  dB  more  signal 
power  than  cooperative  receivers  in  the  binary  case,  while  the  signal  power  difference 
for  BER  =  1.0  •  ICE2  is  ~  6  dB  for  the  4-ary  constellation.  In  both  cases,  the  BER 
gap  between  the  cooperative  and  non-cooperative  receivers  increases  as  Eb/N0  grows. 
In  the  4-ary  case,  the  rate  of  change  of  non-cooperative  BER  improvement  at  high 
Eb/N0  is  near  zero,  suggesting  that  the  ISI  due  to  symbol  dither  is  dominating  the 
error  rate  and  decreasing  the  noise  power  has  little  effect  on  BER.  Because  the  basic 
non-cooperative  receiver  is  outperforming  the  blind  equalization  receivers,  ISI  due  to 
symbol  dither  does  not  appear  to  be  mitigated  in  any  of  the  non-cooperative  receivers 
posited  in  this  thesis.  The  ISI  due  to  symbol  dither,  then,  sets  a  lower  limit  on  the 
non-cooperative  receiver  BER  achievable  for  each  constellation  and  dither  range.  This 
limit  can  be  estimated  empirically  via  simulation  by  removing  the  AWGN  component 
in  the  channel  model.  This  exercise  produces  a  best-case  BER  =  5  •  10-5  for  the  non- 
cooperative  receivers  in  Figure  15.151  (not  shown  on  this  plot)  and  BER  =  4.9  •  10~3 
for  the  non-cooperative  receivers  in  Figure  15.161  This  lower  limit  characterizes  the 
minimum  level  of  security  for  a  particular  set  of  PDW  parameters. 
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Figure  5.15:  Receiver  performance  for  binary  DPSK,  SRRC  (3  =  0.25,  and  dither 
range  0.6 Tsym. 


Figure  5.16:  Receiver  performance  for  4-ary  DPSK,  SRRC  /3  =  0.25,  and  dither 
range  0.6 Tsym. 


5-16 


Figure  5.17:  Receiver  performance  for  binary  DPSK,  SRRC  f3  =  0.25,  and  dither 
range  0.8 Tsym. 

To  degrade  the  non-cooperative  receiver  BER  lower  limit,  the  symbol  dither 
range  can  be  increased.  Setting  the  symbol  dither  range  to  a  higher  fraction  of  the 
symbol  period  also  causes  the  non-cooperative  receivers  to  approach  the  BER  limit 
at  lower  Eb/N0.  In  Figure  [5. 171  non-cooperative  receiver  demodulation  of  the  binary 
DPSK  signal  approaches  a  lower  limit  of  BER  =  4.8  •  10~3  at  about  Eb/N0  =  17  dB 
for  a  symbol  dither  range  of  0.8 Tsym.  The  Eb/N0  required  to  achieve  BER  =  1.0  ■  10-2 
is  9  dB  greater  for  non-cooperative  receivers  than  for  cooperative  receivers.  In  the 

4- ary  constellation,  the  non-cooperative  receiver  lower  limit  is  BER  =  2.4  •  ICE2  and 
is  plotted  in  Figure  5.181 

5- 4  Bandwidth  Considerations 

This  section  compares  the  power  spectral  densities  (PSD)  of  simulated  base¬ 
band  DPSK  waveforms  with  and  without  symbol  timing  dither.  For  standard  DPSK 
modulation,  (without  symbol  dither),  the  width  of  the  PSD  is  dependent  on  the  sym¬ 
bol  period  and  the  pulse  shape.  Figure  157T91  shows  the  effect  of  [3  on  the  null-to-null 
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Figure  5.18:  Receiver  performance  for  4-ary  DPSK,  SRRC  (5  =  0.25,  and  dither 
range  0.8 Tsym. 

bandwidth  using  the  SRRC  pulse  shape.  As  expected,  a  narrow  SRRC  pulse  shape 
in  the  time  domain  corresponds  to  a  wide  null-to-null  bandwidth. 

Because  the  symbol  dither  strategy  does  not  modify  the  pulse  shape  of  the 
transmitted  symbols,  the  spectrum  of  each  symbol  is  no  different  between  a  traditional 
DPSK  waveform  and  the  PDW.  The  time  dither  of  the  PDW  induces  a  phase  shift  in 
the  symbols  in  the  frequency  domain,  but  the  PSD  should  be  consistent  with  standard 
DPSK. 

In  fact,  the  PSD  of  the  dither  waveforms  are  indistinguishable  from  standard 
DPSK  signalling.  In  Figure  5.20[  a  binary  PDW  with  a  dither  range  of  Tsym  (blue 
curve)  produces  a  PSD  similar  to  a  standard  DPSK  waveform  with  dither  range  of 
0  (green  curve  n  Figure  [5.19) .  The  PSD  invariance  with  respect  to  dither  range  was 
observed  for  all  SRRC  pulse  shapes  and  constellations  used  in  this  thesis. 
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Figure  5.19:  PSD  for  standard  binary  DPSK  using  SRRC  pulse  shapes. 


Figure  5.20:  PSD  for  binary  DPSK,  SRRC  pulse  shape  with  (5  =  0.25.  Injecting 
dither  into  the  waveform  does  not  impact  the  average  PSD. 
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VI.  Conclusions  and  Recommendations 


6. 1  Conclusions 

This  thesis  demonstrates  a  communication  system  that  leverages  symbol  timing 
dither  to  hinder  interception  and  exploitation  by  non-cooperative  receivers.  Two 
cooperative  receiver  structures  are  shown  to  reliably  recover  the  digital  information 
with  only  a  slight  bit  error  rate  (BER)  penalty  relative  to  the  optimum  performance 
for  the  DPSK  modulated  waveform.  The  timing  dither  in  DPSK  transmitted  symbols 
causes  ISI  which  is  not  mitigated  in  the  notional  non-cooperative  receivers  simulated 
in  this  thesis.  The  pulse  shape  of  the  DPSK  waveform  and  the  symbol  dither  range 
determine  the  severity  of  ISI  and  set  a  lower  bound  on  BER  of  the  non-cooperative 
receivers.  For  the  signal  parameters  proposed  in  this  thesis,  non-cooperative  receivers 
required  2-8  dB  more  signal  power  than  cooperative  receivers  to  achieve  a  bit  error  rate 
of  1.0%.  For  three  of  the  four  waveform  parameter  sets  analyzed  in  Section  5.3[  non- 
cooperative  receivers  produced  bit  error  rates  in  excess  of  0.1%,  even  for  arbitrarily 
large  values  of  E^/N0. 

Both  cooperative  receivers  use  signal  processing  techniques  to  recover  the  signal 
from  evenly-spaced  samples  of  a  traditional  matched  filter  receiver.  The  modified 
MSE  cooperative  receiver,  which  requires  full  knowledge  of  and  synchronization  to 
the  symbol  timing  dither,  yields  marginally  better  bit  estimation  accuracy  than  the 
LMS  bank  cooperative  receiver.  The  MSE  receiver’s  BER  benefit  is  computationally 
expensive  and  requires  a  matrix  inverse  operation  for  each  symbol  (assuming  the 
symbol  timing  dither  sequence  does  not  repeat).  The  LMS  bank  receiver  assumes  the 
symbol  dither  pattern  repeats,  and  requires  the  period  to  be  short  enough  that  the 
dither  pattern  repeats  many  times  (5,000  times  in  simulations  here)  for  the  equalizers 
to  be  trained.  Thus,  the  low  computational  complexity  of  this  approach  comes  at  the 
cost  of  significant  equalizer  training  periods. 

Of  the  non-cooperative  receivers  considered,  the  basic  sample-spaced  receiver 
with  best-case  synchronization  proved  to  consistently  achieve  better  BER  performance 
than  the  more  complex  non-cooperative  receiver  strategies.  In  practice,  however, 
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non-cooperative  symbol-level  synchronization  may  not  be  easily  achieved  and  the 
oversampled  non-cooperative  receiver  models  may  more  accurately  reflect  expected 
non-cooperative  receiver  performance. 

6.2  Recommendations  for  Further  Study 

In  completing  this  thesis,  several  avenues  for  further  study  have  been  uncovered. 
These  future  research  opportunities  are  concisely  summarized  below  for  easy  reference. 

1.  Many  variations  of  the  dither  waveform  are  possible.  For  example,  the  pseu¬ 
dorandom  dither  applied  to  each  symbol  always  followed  a  uniform  probability 
density  function  in  this  work.  Any  number  of  more  exotic  distributions  could 
be  used.  Also,  simulation  was  restricted  here  to  DPSK  modulation.  The  dither 
strategy  proposed  here  should  be  applicable  to  most  direct  sequence  modula¬ 
tion  schemes,  including  Quadrature  Amplitude  Modulation  (QAM)  and  Pulse 
Amplitude  Modulation  (PAM). 

2.  Simulation  results  were  restricted  to  line-of-sight  (LOS)  channels.  Further  in¬ 
vestigation  could  consider  cooperative  receiver  performance  in  the  presence  of 
multipath  interference. 

3.  The  cooperative  receiver  designs  proposed  here  are  restricted  to  uniform  sam¬ 
pling  of  the  matched  filter  output,  with  the  acknowledgement  that  some  loss  of 
bit  energy  is  expected.  A  more  complex  modified  MSE  architecture  could  use 
complete  knowledge  of  the  dither  sequence  to  sample  the  matched  filter  output 
more  intelligently. 

4.  Lastly,  only  a  few  notional  non-cooperative  receivers  were  used  in  this  thesis. 
Further  work  could  consider  many  other  non-cooperative  receiver  designs.  In 
particular,  the  symbol  dither  estimation  algorithm  could  be  improved.  Also, 
the  impact  of  lion-synchronous  sampling  on  the  non-cooperative  receivers  could 
be  characterized. 
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